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ABSTRACT 

We construct stationary global configurations of both aligned and unaligned logarith- 
mic spiral perturbations in a composite disc system of stellar and isopedically magne- 
tized gaseous singular isothermal discs (SIDs) coupled by gravity. Earlier models are 
generalized to a more general theoretical framework. The thin gaseous SID is threaded 
' across by a vertical magnetic field with a constant ratio of the surface gas mass 

, density to B^. In reference to SID models of Shu & Li, Shu et al., Lou & Shen, Lou 

&: Zou, Shen, Liu & Lou, there exist two classes of stationary magnetohydrodynamic 
(MHD) solutions with in-phase and out-of-phase density perturbations here. Relevant 
parameter regimes are explored numerically. For both aligned and unaligned cases 
Q ' with azimuthal periodicities |m| > 2 (m is an integer), there may be two, one, and no 

^ I solution situations, depending on the chosen parameters. For the transition criteria 

C/3 . from an axisymmetric equilibrium to aligned secular bar-like instabilities, the corre- 

' sponding T/|W| ratio can be much lower than the oft-quoted value of T/\W\ ~ 0.14, 

where T is the total rotational kinetic energy and W is the total gravitational poten- 
tial energy plus the magnetic energy. The T/|>V| ratios for the two sets of solutions 
in different ranges are separated by m/(Am + 4). For the unaligned cases, we study 
marginal stabilities for axisymmetric (m = 0) and non-axisymmetric (m ^ 0) distur- 
bances. By including the gravitational influence of an axisymmetric dark matter halo 
on the background, the case of a composite partial magnetized SID system is also 
examined. The global analytical solutions and their properties are valuable for testing 
and benchmarking numerical MHD codes. For astrophysical applications to large-scale 
galactic dynamics, our model analysis contains more realistic elements and offers use- 
ful insights into the structures and dynamics of disc galaxies consisting of stars and 
magnetized interstellar medium. In particular, in the presence of star burst activities, 
supernovae, hypernovae, superbubbles etc., our open magnetic field geometry in disc 
galaxies bears strong implications on circumnuclear and spiral galactic winds. 

Key words: galaxies: kinematics and dynamics — galaxies: spiral — galaxies: struc- 
ture — ISM: general — MHD — waves. 



1 INTRODUCTION 

The theoretical magnetohydrodynamic (MHD) disc model 
we set out to formulate in this paper is to explore possi- 
ble large-scale global perturbation structures and stationary 
MHD density waves (Fan & Lou 1996; Lou & Fan 1998) in a 
composite system of a stellar disc and an isopedicallly mag- 
netized gaseous disc intended for the interstellar medium 



(ISM). These two gravitationally coupled discs are approx- 
imately treated as 'fluid' and 'magnetofluid' respectively, 
and are both geometricaUy idealized as razor-thin singular 
isothermal discs (SIDs) with the gaseous SID being threaded 
across by an almost vertical magnetic field throughout. In 
astrophysical contexts of large-scale structures in disc galax- 
ies, we also include gravitational effects of a massive axisym- 
metric dark matter halo and adopt a background composite 
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system of two coupled partial SIDs (Syer & Tremaine 1996; 
Shu et al. 2000; Lou 2002; Lou & Shen 2003; Shen & Lou 
2003, 2004a, b; Lou & Zou 2004, 2005; Shen, Liu & Lou 
2005). Our motivation is to construct solutions with com- 
bined analytical and numerical techniques, to understand 
their basic properties, to provide observational diagnostics, 
and to reveal or speculate physical implications. 

ChaJcrabarty, Laughlin & Shu (2003) studied substruc- 
tures in grand-design spiral galaxies, such as branches, spurs 
and feathers, using a two-component disc model in which the 
gas component responds passively and nonlinearly to the 
potential of a rigidly rotating spiral structure involving old 
stars and halos. We here treat a dynamically coupled two- 
component disc system without or with a massive dark mat- 
ter halo and focus on stationary global MHD density wave 
configurations. Specifically, we construct stationary MHD 
configurations for aligned and unaligned logarithmic spiral 
perturbations in a composite disc system of two SIDs with 
flat rotation curves. For observational diagnostics of nearby 
disc galaxies, we also examine phase relationships among 
perturbation patterns of the stellar surface mass density, the 
ISM surface mass density and the isopedic magnetic field. 

We now proceed to provide the more general back- 
ground information relevant to the idealized MHD compos- 
ite disc problem to be formulated and investigated here. 

In a pioneering work on a composite disc system of stel- 
lar and gaseous discs dynamically coupled by gravity, Lin & 
Shu (1966, 1968) proposed a combined approach involving 
a distribution function for the stellar disc and a fluid de- 
scription for the gas disc to derive and analyze the local 
dispersion relation of coplanar galactic spiral density waves 
in the WKBJ approximation. The basic physical scenario 
is that stars form out of gas clouds in the ISM disc, lead- 
ing to the coexistence of a stellar disc and a magnetized 
ISM disc at a later evolution stage of disc galaxy. Since the 
seminal work of Lin & Shu, there have been extensive theo- 
retical researches on density wave oscillations, perturbation 
configurations and stability properties of a rotating compos- 
ite disc system, mainly in the galactic context. Kato (1972) 
studied oscillations and overstabilites of density waves us- 
ing a formalism similar to that of Lin & Shu (1966, 1968). 
Jog & Solomon (1984a. b) discussed the growth of local ax- 
isymmetric perturbations using a two-fiuid formalism in a 
composite disc system. Bertin & Romeo (1988) investigated 
the role of a gas disc for spiral modes in a two-fluid model 
framework. The influence of interstellar gas on oscillations 
and stabilities of spheroidal galaxies was studied by Vander- 
voort (1991a, b). In order to account for the effects of the 
disc thickness, Romeo (1992) adopted a two-fluid approach 
to investigate a two-component disc system with finite disc 
thickness. Lowe et al. (1994) performed an extensive analy- 
sis for morphologies of disc galaxies. Different effective Qef / 
parameters (Safronov 1960; Toomre 1964) have been sug- 
gested for the axisymmetric stability of a composite disc 
system in a two-fluid formalism by Elmegreen (1995) and 
Jog (1996). Lou & Fan (1998b) recently used the two-fluid 
formalism to study properties of open and tight- winding spi- 



ral density-wave modes in a composite disc system. Lou & 
Shen (2003) discussed stationary global perturbation struc- 
tures in a two-fluid formalism and offered a more straightfor- 
ward D— criterion for the axisymmetric instability in a com- 
posite SID system instead of a redefinition of a new Qeff 
parameter (Shen &: Lou 2003). Considering the magnetic 
field embedded in the gaseous disc, Lou & Zou (2004, 2005) 
studied the stationary global coplanar MHD perturbation 
structures and axisymmetric stability in a composite SID 
system. 

A series of astrophysical disc problems involves stability 

analysis of a SID. Since the pioneer work of Mestel (1963) 
more than four decades ago, numerous theoretical and nu- 
merical studies have been carried out in this subject area 
(e.g. Zang 1976; Toomre 1977; Lemos, Kalnajs & Lynden- 
Bell 1991; Lynden-Bell & Lemos 1999; Goodman & Evans 
1999; Chakrabarti, Laughlin & Shu 2003). An important 
breakthrough by Syer & Tremaine (1996) was to derive the 
semi-analytic solutions for stationary coplanar perturbation 
configurations in a class of power-law discs (i.e., SID is only 
a special case). Shu et al. (2000) obtained stationary solu- 
tions for global perturbation configurations in an isopedi- 
cally magnetized SID with a flat rotation curve. Based on 
extensive numerical calculations, they interpreted these sta- 
tionary aligned and unaligned logarithmic spiral conflgura- 
tions as the onsets of bar-type and barred-spiral instabilities 
(see also Galli et al. 2001). As a different yet complementary 
work to the analysis of Shu et al. (2000), Lou (2002) stud- 
ied a coplanar MHD perturbation analysis in a single back- 
ground SID embedded with an azimuthal magnetic field, 
from the perspective of stationary fast and slow MHD den- 
sity waves (FMDWs and SMDWs; Fan & Lou 1996; Lou & 
Fan 1998a). Lou (2002) also derived a form of magnetic virial 
theorem for an MSID and suggested by analogy that the ra- 
tio of the disc rotational kinetic energy to the sum of the 
gravitational and magnetic energies may be crucial for the 
onset of bar-like instability in an MSID system. To model 
large-scale effects of a galactic magnetic field, it would be 
more realistic to treat large-scale structures and dynamics as 
a composite disc system of stars and a magnetized ISM. Lou 
& Shen (2003) made an foreray into this model problem, con- 
structed stationary aligned and unaligned logarithmic spiral 
configurations in such a composite SID system and further 
examined axisymmetric instability properties (Shen & Lou 
2003). Meanwhile, Lou & Zou (2004, 2005) incorporated a 
coplanar magnetic field into the composite SID system for 
a more specific analysis. 

The geometry of an isopedic magnetic field proposed 
by Shu & Li (1997) is an open magnetic field configura- 
tion across the disc plane. For disc galaxies, open magnetic 
fields are most likely interlaced or intermingled with 'closed' 
coplanar magnetic fields in a stochastic manner.* For a much 
simplified theoretical model analysis at this stage, we take 



* The closed and open magnetic fields over the solar surfax;e into 
the lower solar corona serve as an empirical analogy. 
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an isopedic magnetic field geometry all over a composite 
SID system. Together with other physical mechanisms (e.g., 
massive star formation activities, supernovae, hypernovae, 
hot bubbles and superbubbles etc.), such an open magnetic 
field geometry can generate and channel hot gas outflows, of- 
ten referred to as 'galactic winds' under increasing scrutiny 
of multi-band observations and the X-ray band in particu- 
lar. Along with observational and theoretical studies of stel- 
lar winds (e.g.. Burke 1968; Holzer & Axford 1970), an ini- 
tial impetus was made to formulate the seminal concept of 
galactic-scale winds blowing away from both sides of the 
disc plane (e.g., Johnson & Axford 1971). As galactic winds 
are expected to be a relatively weak phenomenon in terms 
of radiative signatures, we start to accumulate some sub- 
stantial observational evidence in recent years. It is gener- 
ally conceived that supernova explosions and stellar radia- 
tions drive galactic outflows. Strickland, Ponman & Stevens 
(1997) have reported X-ray observations of galactic outflows 
from the starburst galaxy M82, which may imply an over- 
all configuration of a galactic wind. On both observational 
and theoretical grounds, magnetic field configurations away 
from a galactic plane should play dynamical, diagnostic and 
geometric roles in probing galactic winds. As a result of such 
open magnetic field, the outflow of hot gas can lead to large- 
scale observable effects. 

For such a magnetized composite disc system, MHD will 
play an indispensable role in the gaseous MSID and reveal 
more realistic aspects of coupled large-scale dynamics as well 
as diagnostic information. These important MHD disc prob- 
lems (Shu et al. 2000; Lou & Shen 2003; Lou & Zou 2004, 
2005; Shen, Liu & Lou 2005) are not only interesting for 
their own sakes, but also serve as necessary steps for devel- 
oping more and more realistic physical models. We construct 
here stationary perturbation configurations for aligned and 
unaligned logarithmic spiral cases in a composite system of 
a stellar SID and an isopedically magnetized gaseous SID, 
and discuss their stability properties. For numerical MHD 
simulations into the nonlinear regime, the global perturba- 
tion solutions here are valuable for initializing, benchmark- 
ing and understanding the MHD code development. 

We adopt a relatively simple formalism for a compos- 
ite disc system containing infinitely thin fluid (for stars) 
and magnetofluid (for the magnetized ISM) discs dynam- 
ically coupled by gravity. We describe in Section 2 the iso- 
pedic magnetic fleld and a few important theorems for the 
gaseous MSID. In Section 3, we determine conditions for the 
background axisymmetric equilibrium state for both stel- 
lar SID and gaseous MSID, and derive linearized equations 
for MHD perturbations. The global MHD perturbation so- 
lutions can be classifled as aligned perturbations and un- 
aligned logarithmic spiral perturbations; these solutions are 
analyzed in detail in Section 4. The exact solutions of sta- 
tionary MHD perturbations, their properties and their corre- 
sponding phase relationships among perturbation variables 
are examined and summarized in Section 5. Several details 
of mathematical analysis are included in Appendices A and 
B for the convenient of references. 



2 ISOPEDIC MAGNETIC FIELDS 

Magnetic fields may be generated and sustained through 
nonlinear MHD dynamo processes in an electrically con- 
ducting gas medium in motion (Parker 1979 and extensive 
references therein). Several numerical simulation results and 
some theories for galaxy formation lend supports to the ba- 
sic assumption that the mass-to-magnetic flux ratio A may 
be constant in space and time. For example, model calcula- 
tions and numerical simulations of cloud core formation by 
ambipolar diffusion tend to produce gravitationally unstable 
inner regions with approximately constant values of A (e.g., 
Nakano 1979; Lizano & Shu 1989). In the models of Basu & 
Mouschovias (1994), the mass-to-magnetic flux ratio A re- 
mains almost constant when the variation ranges of density 
and magnetic field fiux are large. In the theoretical formu- 
lation of Shu & Li (1997) and Shu et al. (2000), it was thus 
presumed that A is a constant; in particular, they referred 
to these earlier results to justify this assumption. Here, we 
emphasize the theoretical fact that A remains constant is a 
natural consequence of the standard ideal MHD equations. 
In cylindrical coordinates (r, ip, z) and by using the gas mass 
conservation equation 12011 below and the £— component of 
the magnetic induction equation 
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in the thin disc geometry, one can readily show that 
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dt\Bj ' dr\B,J ' r dp\B, 
where the relevant notations bear the conventional meanings 
in MHD. In reference to the works of Mouschovias (1994) 
and Li & Shu (1996) as well as the earlier numerical results 
(Nakano 1979; Lizano & Shu 1989), it is fairly clear that the 
above equation provides a natural explanation for the fact 
that A remains constant in the evolution of an MHD disc. 

For the model problem under consideration, we ignore 
vertical structures within the disc thickness as the horizontal 
scale of the problem is very much larger. To analyze proper- 
ties of a drastically flattened disc, we adopt the cylindrical 
coordinates (r, ip, z) in our mathematical model prescription. 
We assume that the mass density p and the electric current 
density je in our model to have nonzero values only within 
a narrow range Az <C r about the disc midplane z = Q. We 
are justified to define a surface mass density E and a surface 
electric current density J by 



and 



T.{r,p,t) = 



3{r,p,t) 



/ p{r,p,z,t)dz 

J Az 



je{r,p,z,t)dz 



with J — JrGr + J^e^p having only two orthogonal com- 
ponents coplanar with the disc and the component of J 
within the disc thickness is omitted. 

As for the two SIDs, one is a stellar SID approximated 
as a fluid for large-scale dynamics and the other is a magne- 
tized gas SID treated as a magnetofluid. Over large scales, 
we ignore the interaction between the magnetic field em- 
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bedded in the ISM and the stellar disc. It then comes the 
physical scenario that the stellar and magnetized gas discs 
are dynamically coupled by gravity, while the magnetic field 
directly interacts with the gaseous ISM disc. 

On the basis of the Poisson equation relating the gravi- 
tational potential and the mass density, there will also be a 
jump in the vertical gravitational field given by 47rGE across 
a thin disc. In other words, the vertical gravitational field 
just above and below the disc point towards the midplane 
z = and have the same strength of 27rGE. We choose to 
define the ratio t) of the horizontal gravitational acceleration 
/ll (continuous across the thin disc in vertical direction) to 
the vertical gravitational acceleration just above the disc as 
a dimensionless parameter, namely 

' (27rGE) 

and the local value of the mass-to-magnetic flux ratio as 

S"(r-,y,t) 
B,(r,ip,Q,t) ■ 

We then take A to be constant in both space and time, 
when computed at the footpoint of a magnetic field line 
anchored through the disc (the adjective 'isopedic' comes 
from the Latin word 'ped' for foot). For ideal MHD with the 
above frozen-in condition, an initially isopedic disc retains 
the same value of A throughout the entire MHD evolution. 

It is convenient and advantageous to nondimensionalize 
this ratio in the form of 

A = 27rG*A . 

We now discuss the two MHD theorems proven by Shu & Li 
(1997) that form the fundamental part of the MHD equa- 
tions to come in our theoretical model analysis. 

THEOREM 1: For an arbitrary distribution of the sur- 
face mass density Jl{r,ip,t), which may or may not be in 
mechanical equilibrium in the horizontal directions, the mag- 
netic tension force in the plane of a thin isopedic disc equals 
— times the horizontal self- gravitational force /y . 

THEOREM 2: // we also assume that the disc remains 
in a vertical magnetostatic equilibrium, then the magnetic 
pressure integrated over the disc thickness may be approxi- 
mated as (1 -I- ri'^)/{X^ + r]'^) times the thermal gas pressure 
integrated over z. 

By the first theorem, the magnetic tension force /ten 
tangential to the disc plane is simply 

/ten = -f\\l'>^ 

and the combined forces of magnetic tension and gravita- 
tional field act in a form of a diluted horizontal gravity 

/ = /ten + f\\ = e/|| , 

where e = 1 — is the reduction factor. 

By the second theorem, it follows approximately that 

where Hm is the projected magnetic pressure acting on the 
gas material and Hg is the vertically integrated thermal gas 



pressure. Therefore, the total projected pressure is 

H = H„ Hg = OHg = Qa^Yf , 

where the enhancement factor O is given by 

(1 + A^+V) 
" (A2+r72) • 

Coming to a composite disc system of two gravitationally 
coupled discs with the gas SID being isopedically magne- 
tized, the above two theorems should be modified accord- 
ingly. By carefully examining the two theorems in details, 
we note that as the stellar disc does not interact with the 
magnetic field, the first theorem of Shu & Li should retain 
its basic form with some relevant parameters being prop- 
erly adjusted. For example, A should now stand for E*/_Bz 
where is the surface mass density of the gas disc and 
the magnetic tension force /ten acts on the magnetized gas 
SID. Naturally, we introduce an appropriate r\ parameter 
for convenience, namely r\ = /||g/(27rGE®), where /y^ is the 
contribution to the total /y from the gaseous MSID. For the 
two SIDs having the same structure and density profile, we 
further obtain 

- f\\9 ^ if\\g + f\\s) ^ f\\ 
' 27rGE9 27rG(E9-HE«) 27rGE ' ^' 

where E° is the surface mass density of the stellar SID, E is 
the total surface mass density and /y^ is the contribution to 
the total /y from the stellar SID. 

For the second theorem of Shu & Li, the generalization 
and extension are somewhat different. We now examine the 
procedure of Shu & Li to generalize the second theorem 
for a composite disc system with one stellar SID and one 
isopedically magnetized SID. 

The force per unit volume associated with the magnetic 
pressure gradient is given by -V[B'^/{Stv)], where 

= B^r +Bl + Bl . (2) 

By integrating — V[-B^/(87r)] from z — — oo to z — -\-oq, the 
contribution from the z— component of the magnetic pres- 
sure gradient vanishes for — > as |z| ^ -|-oo. We thus 
obtain the parallel magnetic force per unit projected surface 
area as 

-Vy / j—dz , (3) 

where the parallel gradient operation within the disc plane 
is deflned by 

V = e'^ — + ^ — 
" 9r r dtp 

Away from the disc, the total integral of equation iPJ 
does not act on materials contained in the gaseous MSID. In 
the vacuum regions above and below the MSID, the horizon- 
tal force of magnetic pressure gradient does not vanish but 
counteracts instead a nonzero force of magnetic tension, such 
that the sum produces a force-free environment in vacuum 
(i.e., no electric current flowing around). We specifically de- 
note the part of the projected magnetic pressure acting on 
the gaesous MSID as Ilm- For a magnetized gas disc of a 
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thin but nonvanishing thickness, we estimate Hm by 

= / —dz = , (4 

Jaz Stt (47r) 

where zq is the effective half-thickness of the gaseous MSID 
and -6+ is the value of evaluated just above the upper 
disc surface. While our estimate for Ilm is plausible on the 
dimensional ground, the possibility remains for a multiplica- 
tive factor of order unity whose exact value would depend 
on the vertical structure within a thin MSID. 

In order to determine the two expressions of and 
we begin with a generalization which states that a quasi 
vertical magnetostatic equilibrium in a thin gas disc requires 
the sum of the gas, magnetic and gravitational pressures to 
be independent of the vertical height z, namely. 



(5) 



where Pg is the thermal pressure of the gaseous MSID and K 
is an integration constant; and by assuming that the gaseous 
MSID and the stellar SID interact only via the mutual grav- 
ity, we have the total effective surface mass density in equi- 
librium as 



(Ps + P3)dz 



(6) 



In the midplane 2: = 0, we have by symmetry cr = 0, Br = 
B^ = 0, Bz = E3/A, while Pg = ngl{2zo) with Hg being 
the vertically integrated thermal gas pressure and zo is the 
half-scaleheight in the vertical direction. Thus equation 
yields at 2; = the following expression for constant K, 



9\2 



(7) 



220 87rA2 

Meanwhile for r ^ 2: ^ zq, we have Pg = 0, Ps = 0, a = 
E" -I- E", and B = {T,o/A)e^ - f||g/(27rG'A) [see derivations 
of Shu & Li (1997) that lead to their equations (2.14) and 
(2.15)]. In these spatial regions, equation @ leads to 



87rA2 327r3G2A 



2 + Jg(e«' + e^)' 



(8) 



By introducing the gas-to-stellar surface mass density ratio 
S = E^/E" from the two SIDs, we derive from equations Q 
and ||HJ| the following relation 



Hi 

zo 



7rG(E 



9\2 



+ 



(47rGA2) 



g,,[il + Sr)^±frf] 



It then follows immediately that 



20 = 



7rG(E9)2 [(1 -f 5)2A2 + 52j,2] 

Just above the disc, equation leads to 



Bl 



K 



-G(E 



9\2 



1 



2 -V- > 

Using expression Q for constant K, we then obtain 



Bi 
Stt 
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Yzo 



g , ^G(E 
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7rG(E' 
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2A2 



'4 



(9) 



(10) 



(11) 



(12) 



Multiplying equation 11211 by 22o with 20 being given explic- 
itly by expression 1101 . we derive from definition Q 

2oB2 _ {l + rf)^g 



(47r) [(H-5-i)2A2 +772] • 

It is convenient to introduce a modified A parameter A 
(1 -I- (5~^)A such that expression 1131 becomes 

(1 + '?') 



(13) 



n. 



Ho 



(14) 

(A2+r;2) 

By summing up the gas and magnetic pressure contribu- 
tions, we obtain the total pressure in the form of 

n = + n„ = QUg , (15) 

where the modified enhancement factor Q is given by 

(l + P + 277^) 



e 



(A2 



(16) 



This completes our extension and generalization of the sec- 
ond theorem of Shu & Li as applied to the MHD problem 
of a composite disc system consisting of one stellar SID and 
one gaseous isopedic MSID. 

At the end of this section, we summarize the two ex- 
tended theorems for a composite MSID system with an iso- 
pedic magnetic field as follows. In the presence of an isopedic 
magnetic field, the horizontal self-gravity force f||g and the 
gas pressure integral Iig{= OgE^) of the gaseous MSID will 
be modified as 



Ha 



Hg 



where e = 1 — A ^ and is defined above by equation 1161 
that involves essentially three dimensionless parameters 5 = 
EVE^ A = 2^Gi/2EVB,, = /||g/(27rGE«) = /||/(27rGE) 
and A = {l + 5-^)X. 



3 FLUID-MAGNETOFLUID FORMALISM 
FOR TWO THIN COUPLED SIDS 

3.1 Basic Equations for a Composite MSID 
System and the Background Equilibrium 

The two coupled SIDs located at 2 = are both approx- 
imated as infinitesimally thin for expediency and they are 
coupled through the mutual gravity. For large-scale station- 
ary MHD perturbations, diffusive processes such as viscosity, 
ambipolar diffusion and thermal diffusion etc. are all ignored 
in our formulation for simplicity. For physical variables un- 
der consideration, we use either superscript or subscript s 
to indicate the association with the stellar SID and either 
superscript or subscript g to indicate the association with 
the gaseous MSID. In the cylindrical coordinates (r, ip, 2), 
the ideal fully nonlinear fluid-magnetofluid equations for a 
composite MSID system can be readily written out. 

For the stellar SID in the 'fluid' approximation, we have 



aE" 



1 3 , ^ s 
r or 



r2 d(p 



(17) 



6 Y.-Q. Lou & Y. Wu 



for the two-dimensional stellar mass conservation, 



■ + u — h 



dt 



dr 



dip 



(18) 



for the radial momentum equation, and 

df ^sdf df 1 d 2 
-qT^ — + Z2 — = — («»^ ) 



dr 



T," dip' 



dip 



(19) 



for the azimuthal momentum equation, where the 'isother- 
mal' approximation has been invoked. 

In parallel, for the gaseous isopedic MSID in the mag- 
netofiuid approximation, we have 

+ -^(^SV) + 47r^(SV) =0 (20) 



dt 



r dr 



dip 



for the two-dimensional gas mass conservation. 



du 



J du^ du^ j 



92 



1 d 



dt 



dr 



dp 



E9 dr' 



for the radial momentum equation, and 
df _^^sdj^_^fdf _ 1 d 



dt 



dr 



dip 



d(p 



(ea^E" 



a(0^+e^ 
' dr 

(21) 

dp 



(22) 

for the azimuthal momentum equation, where e = 1 — 1/A^ 
and the modified here is defined by equation 11611 . The 
dynamical coupling between the two sets of fiuid and mag- 
netofluid equations is caused by the gravitational potential 
through the Poisson integral relations 

-GJ:'{r',Tp,tydr' 



(j>"{r,ip,t) 

and 

<l>'(r,^,t) 



J.2 _ 2rr' cosixp — p)Y^'^ 



(23) 



-GY.o{r',ilj,t)r'dr' 
]j.i2 ^j.2_ 2rr' cos(i/) — ip)Y^'^ 



(24) 

as well as through the effects of 5 and r) in the Ilm term. In 
equations 1171 — 1241 1. E* is the stellar surface mass density, 
u" is the radial component of the bulk stellar 'fiuid' velocity, 
is the z— component of the specific angular momentum 
of the stellar 'fiuid', and as is the stellar velocity dispersion 
presumed to be constant (or an equivalent 'isothermal sound 
speed' for the stellar SID), OgE* stands for an effective pres- 
sure in the isothermal approximation. For physical variables 
of the gaseous isopedic MSID, we simply replace the super- 
script s by g systematically. Here, we assume that the stellar 
and magnetized gaseous SIDs interact mainly through the 
mutual gravitational coupling on large scales. As the isope- 
dic magnetic field only interacts with the gas disc, and e 
only appear in the MHD equations for the gaseous MSID. 

Before an MHD perturbation analysis, one needs to 
establish a background magneto-rotational equilibrium for 
the composite MSID system that is dynamically consistent 
with equations 11711 — 12411 . For that purpose, we presume 
an axisymmetric background for both stellar and magne- 
tized gaseous SIDs, with the same form of power-law surface 
mass densities (i.e., E oc r~^) yet with different fiat rotation 
curves in general. 

Using the radial momentum equations 11811 and I2H for 
the background equilibrium of axisymmetry with Uq = = 



Q,Vls= jo/r"^, and VI g 



and 



r 



-Q.gr = e- 



jQ/r'^, we readily obtain 
27rG(Eg + Eg) 



27rG(Eo + eY.1) 



(25) 



(26) 



where 5 = Eq/Eq, and Q,s = UsDs/r and 0,g = UgDg /r 
for the two rotational Mach numbers Ds and Dg, respec- 
tively. The two equilibrium surface mass densities Eq and 
Eq can be immediately determined as 



E^ = 



and 



ya 



2-KGr (1 + 5) 



Qal 5{l + Dl 



(27) 



(28) 



2-KGr {l + eS) ' 
For stationary global MHD perturbation configurations to 
persist, Ds and Dg can only take on specific values (e.g. Lou 
& Shen 2003; Shen & Lou 2003, 2004a, b; Lou & Zou 2004, 
2005). In the above equations, the two parameters 5 and e 
are constant, while O is a variable in general. However, in 
our MHD perturbation analysis, O variation will be ignored 
merely for simplicity without losing essential physics (see 
Shu & Li 1997; Shu et al. 2000). We now further introduce 
a new dimensionless parameter (3 = a^sl'^g for the square of 
the 'sound speed ratio'. 

Because |/||g| = |d(^Q/dr| = 2-kGYPq in an axisymmetric 
equilibrium state, we take r] = \ and have O > 1 in gen- 
eral. For a positive Eq as a necessary physical requirement, 
the equilibrium solution 128t implies the following three in- 
equalities in sequence 



H-e5=l+|l-^)<5>0 



5 + l>^, 

A(l + <5-i)>i . 

As 5 > by definition, it is obvious that \(&~^ -I- 1) > A. 
A multiplication of this inequality with the third inequality 
above on both sides leads to A = A(l + 5^^) > 1 for a physi- 
cal magneto-rotational background equilibrium. It then fol- 
lows that the enhancement factor defined by equation 1161 
should be constrained within a finite range of 1 < O < 2. 

We now discuss the parameter e. In Shu & Li (1997), 
they require e > which is the same as A > 1 for a gravita- 
tionally bound disc. In the case of a single SID, the physical 
explanation of this condition is that a rotating disc with 
self-gravity must be held together in a radial force balance. 
As a result of rotation and singular isothermal density dis- 
tribution, the centrifugal force and the effective and thermal 
pressure forces all point radially outward. The only counter- 
force to keep the disc in equilibrium is the radially inward 
gravity force. For e < 0, the magnetic tension force will be 
stronger than the horizontal gravitational force; in this case, 
there will be no rotational equilibrium for the disc system 
at all to begin with. In our composite model containing one 
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SID and one isopedic MSID, the necessary condition for an 
equilibrium to exist is 1 + 6(5 > in comparison with the 
condition of e > of Shu & Li. More precisely, we should 
have < e < 1 for a single SID and < (1 + e(5) < (1 + 5) 
for a composite system of two gravity-coupled SIDs. 



3.2 MHD Perturbation Equations 

We now introduce nonaxisymmetric MHD perturbations to 
the disc system, denoted by physical variables with subscript 
1, in both stellar SID and gaseous isopedic MSID, namely 

E'' = + E'' = E« + E? , 



U = + = lii 



9 9 \ 9 9 



J = Jo+J\ , J - Jo +31 ■ 

It is then straightforward to write down the linearized MHD 
perturbation equations in the forms of 

aE? 1 a 9E? Eg dji 



f +rf^X + ^^.f = -^(a?| + 0? + 0n (31) 
for the stellar SID, and 



aE? , 1 9,^^,,,,, , ^ 9E£^E|ajf 



(32) 



d 



—±+ng-±^2Q/-^ = -— Qa^^+^l+e^f , (33) 



at 



dr 



ya 



2A+rn,ui + n,^ = -4-(eal$+ri + e^l] (34) 



for the gaseous isopedic MSID, with the total gravitational 
potential perturbation = (jil + (f)\ given by the sum of 

,\ L< r -G^\{r' ,^ ,t)r'dr' 
Mr,V,t) = fd^J^ [,,2+,2,2rr'cos(V>-y)]V^ ^^^^ 

-GE?(r' ,ip ,t)r'dr' 



and 



g [r'2 +r2 - 2rr'cos(-i/) - (^)]i/2 ' (3^) 

Assuming a Fourier harmonics in the periodic form of 
exp[i(a;t — mip)] for MHD perturbation solutions in general 
(after taking the real part), we write for perturbations 

(37) 



E'i = {"t) exp[j(Ljt — Tii^j)] 



j[ = J^ir) exp[i(cijt — rmp)] 



= V'{r) exp[i{Lut — rriLp)] , 



(38) 



(39) 



(40) 



where I — s or g denotes associations with stellar or gaseous 
discs respectively. By substituting expressions (1371 — MOI I 
into equations 1291 — II36I I. we derive for the stellar SID 



Id J" 

i{u - mQs)l^" H r(rEgf/") - imEo^ = 

r dr 



for the gaseous isopedic MSID 
r dr 



mn'')r + rflsU" =im{ a'i^ + V' + V' 



J9 

imE« — 



i{iu-mn^)U''-2ng— = --^[eal^+V' + eV 



JO 



d 



2 M 



and for the two gravitational potential perturbations 

—G^i"{r') cos {m^p)r'dr' 
+ — 2rr' cos i/)]-'-''2 ' 

—G^P{r') cos {mtp)r'dr' 



Vir) ^ (b dtp 



V^ir) 



dtp 



+ r2 



2rr' cos i/jJI/^ 



(41) 
(42) 
(43) 

(44) 
(45) 
(46) 

(47) 
(48) 



We use equations 14211 and 1431 to express U" and J* in 
terms of = (a^/^VEg) + V" + for the stellar SID 
and similarly, use equations l|45^ and H46|l to express [7® 
and JO in terms of = (ea^/^f/E^) + V + eV^ for the 
gaseous MSID. In the following analysis, we shall assume for 
simplicity that O remains constant in the presence of MHD 
perturbations in a composite MSID system (see Shu et al. 
2000). The resulting expressions then become 



—i{uj — mils) 
{lj - mQ.sf - 2^1 
and 

J" (ijO - mr2s) 

for the stellar SID, and 



2Qsm 



r{uj — mQs) 
m Q.S 



d_ 

dr 

d 



r (uj — mD,s) dr 



(49) 



(50) 



and 

r {u! — mQg)^ 



-2^2 



20,gm 



2ni 



r{uj — niQg) 
n Q,„ 



d_ 

dr 

d 



^0 (51) 



(52) 



(uj — mQ.g) dr 
for the gaseous isopedic MSID, respectively. 

A substitution of expressions H49|l and H5()|l into equa- 
tion Ij41|l leads to 



1 d { rT,Q{Lo - mQs) 



r dr \ (w - mf7,)2 - 2r22 
mT,(){ij — mO,s) 



2Qsm 



d_ 

dr 



r{u> — mQ,s) 
m Qs d 

r {u! — mQs) dr 



(53) 



r[{u - mQs)^ - 201] 
for the stellar SID. Likewise, a substitution of expressions 
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15111 and 1521 1 into equation 14411 leads to 

1 d ( rEg(tj — mf2g) 
m'Eg{u! — m.Qg) 



20,gm 



d_ 

dr 



r[{u) ~ mQg)^ 



r[uj — mQg) 
m Qg 
r — milg) dr 



d 



5(9 



(54) 



for the gaseous isopedic MSID. Equations 1531 and 15411 
should be solved together with Poisson integral relations 
lITTIl and M . 

For stationary MHD perturbation solutions (u) — 0) 
with zero pattern speed, we readily have 



r dr\ (m2 - 2)Q.s 



r(m2 - 2)Qs 
for the stellar SID, and 



2 + 4 
dr 

'■m? d 
r dr 



(55) 



1 d 



r dr\ {m? - 2)^9 



2 + 4 
dr 



mEg 



m 
r 



d_ 

dr 



(56) 



r{m? - 2)Q.g 

for the gaseous isopedic MSID, respectively. 

Using conditions and relationships 12511 — 12811 . we fur 
ther reduce the above two equations to the forms of 



and 



{m? - 2)Q.lr \ r 



(m2 - 2)Q?gr \ r 



2— r— I*' 

dr dr^ 



dr 



: (57) 



(58) 



for the stellar SID and the gaseous isopedic MSID, respec- 
tively. For nonzero and fig, we have 



1 



Di(m2 - 
X ^r^" + 

for the stellar SID and 
1 



d 



2)\ r dr 
lJ_D^]rj-V^ 
2-kG 1 + (5 



d^ 

dr^ 



(59) 



m 

2)1 r 



r/i^ + 



2 „ r ~dr 



dr'^ 



2.G (1 + .5) JJ " 

for the gaseous isopedic MSID, respectively. Equations 15911 
and are to be solved simultaneously with Poisson inte- 
grals H47|l and l|48|l in order to determine the two rotational 
Mach numbers Di and Di. 



By careful and independent measurements of galactic rota- 
tions (e.g., Rubin 1987; Kent 1986, 1987, 1988), the more 
or less flat rotation curves of galaxies determined by dif- 
ferent methods clearly imply the presence of dark matter 
halos in spiral galaxies. To make our formulation more gen- 
eral and realistic, we here propose a theoretical model of a 
composite partial MSID system to take into account of the 
Newtonian gravitational effect of a presumed axisymmetric 
dark matter halo (e.g., Binney & Tremaine 1987; Berlin & 
Lin 1996; Syer & Tremaine 1996; Shu et al. 2000; Lou 2002; 
Lou & Shen 2003; Shen & Lou 2003, 2004a, b; Lou & Zou 
2004, 2005; Shen, Liu & Lou 2005). For a composite partial 
MSID system, the preceding mathematical formulation for a 
composite MSID system remains to be its essential part yet 
with an important distinctive feature of adding an axisym- 
metric gravitational potential in the background magneto- 
rotational equilibrium. In numerical simulations designed for 
galaxy formation, velocity dispersions of dark matter 'par- 
ticles' in the halo are typically high (e.g., Hoeft et al. 2004). 
For simplicity, we may therefore ignore dynamical pertur- 
bation responses of this massive dark matter halo to MHD 
perturbations in the composite partial MSID system. More 
specifically, we introduce an additive gravity term —d^/dr 
to the basic equations 1181 . 1191 1 and 1211 1. 1221 1. where <I> 
represents the gravitational potential contribution from the 
background dark matter halo of axisymmetry. In the two 
components of the momentum equation for the stellar disc, 
we then have 



du" 
dt 



+ u' 



+ 



dr 
E" dr 



r^ dip 
(a^E^) 



9(<^^ 



+ (/>» + $) 



dr 



(61) 



dt dr r^ dip 



1 d 2 9(0" + 0« + $) 

'^d^'^""^ > d^ 



(62) 



meanwhile, for the two components of the momentum equa- 
tions in the gaseous isopedic MSID, we have 



dr r^ dip 



dt 



+ u- 



E9 dr' 



— ^- u 

dt dr 



df f df 



+ 



dip 



1 d 



djcj}' + e(f>^ + $) 
dr 



djcj)" + e4,<' + $) 
dip 



(63) 



(64) 



3.3 A Composite System of Partial (M)SIDs 

Since the early pioneer work of Zwicky (1933) and Smith 
(1936), evidence has gradually emerged for the presence of 
dark matters in clusters of galaxies. In the 1970s, system- 
atic observations of field spiral galaxies further revealed the 
generic presence of dark matters in isolated spiral galaxies. 



For such a composite partial MSID system, we now in- 
troduce a dimensionless ratio T parameter defined by T= 
{(l)^ + <f)")/{<l)^ + (t>' + ^) (Lou & Shen 2003; Shen & Lou 2003; 
Lou & Zou 2004, 2005). The background magneto-rotational 
equilibrium should be modified accordingly. As before, we 
still write fis = asDs/r, fig = agDg/r, and uq = Uq = Q 
for the background. The background equilibrium equations 
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then become 



2nG 



:sg + Eg) 



and 



r 



2-kG 



(65) 



(66) 



It follows immediately that the two modified equilibrium 
surface mass density profiles are now given by 



and 



ys — 



27rGr (l + S) 
; (1 



(67) 



(68) 



2nGr 1 + [1 - J^(l - e)]S 
respectively. From these two radial equilibrium conditions, 
we readily derive the requirement 

a^{l + [l-.^(l-e)]<5}(l + D,') = ea^(l + 5)(l + Z)^) (69) 

resulting from the same total gravitational potential that 
couples the two SIDs. We refer to such a system as a com- 
posite system of partial (M)SIDs for < < 1, in contrast 
to a composite system of full (M)SIDs with JT = 1. 

The linearized MHD perturbation equations in a com- 
posite system of partial (M)SIDs should take the same forms 
as those in a composite system of full (M)SIDs explicitly 
written down earlier in this section. For stationary perturba- 
tion configurations in the modified equilibrium with u — 0, 
we then have for the stellar partial SID 

/nf_ _^d_ ^ (f 

2)\r dr dr^ 



1 



2nG (1 + &) 
and simultaneously for the gaseous partial MSID 



= 



(70) 



Dl{m? - 2) 



— -2 — 

r dr 



d^ 
dr^ 



2-kG l-f [1 - JP'(1 - e)]5 
where T factor appears explicitly in both relations. 

3.4 Two Different Limiting Procedures 

Here, we discuss specifically the axisymmetric case of m = 
and start from the very beginning with w 7^ in Fourier 
harmonics I37II — 14UII . By substituting expressions 1371 — I4UII 
with m = into equations I29II ~ 13611 . we derive 



r dr 



iujU" 



2Qs— = 
r 



iujr +rQ,sU'' = 



for the stellar SID, 



r dr 



(72) 
(73) 
(74) 

(75) 



JO 



dr 



2n„— = Qal^ +V' + eV 



s 



iujjo + rngUO = 
for the gaseous isopedic MSID, and 

-Gfi"{r')r'dr' 



Vir) ^ f dtp 



Vir) 



dip 



r'2 + r-^ — 2rr' cos 'i/)]-'-/^ ' 



-Gfi^ {r')r' dr' 



(76) 



(77) 



(78) 



(79) 



/o [r'2 -f — 2rr' cos t/iji/^ 

for the two gravitational potential perturbations. 

Again, we use equations 17311 and 1741 to express U" and 

r in terms of = (a^^VEg) + + for the stellar SID 

and in parallel, we use equations l)76|l and H77|l to express 

and in terms of = Oialfi^ /Eg) + + eV^ for the 

gaseous MSID. The resulting expressions are 

iij d^" , , 

U = , . — . (80) 



3' 



(lj2 - 2^1) dr 



r (aj2 _ 2ni) dr 
for the stellar SID, and 



r 



(a;2 - 2^2) dr 



(a;2 -2^2) dr 



(81) 

(82) 
(83) 



for the gaseous MSID, respectively. 

A substitution of expressions 18UII and 18111 into equa- 
tion 1721 1 leads to 

wrEg d* 



= 0;^ -I — 

r dr 



(w2 -2^2) dr 



(84) 



for the stellar SID. Likewise, a substitution of expressions 
H82|l and H83II into equation H75|l leads to 



= , (71) = Ljfi' + 



ld_ 

r dr 



d*9 



-2J72) 



dr 



(85) 



for the gaseous isopedic MSID. Ordinary differential equa- 
tions (ODEs) (IHU and ^ are to be solved with Poisson 
integrals 1781 and 1791 1 simultaneously. 

The case of m = is special. If we set lj — 0, all MHD 
perturbation equations are satisfied in a trivial manner. The 
harmonic factor exp [i{ijt — mip)] becomes unity for both 
ui — and m = 0. In order to construct stationary con- 
figurations, we set a; ^ to derive limiting solutions. 

The two equations above can be cast into 



°^'"r'+(c.2_2n2) 

for stellar disc and 



d 



(a;2 - 2^2 )rdr 



d^ 

dr^ 



= uj{iJ,o + 



ys 



d d2 



(aj2 



2ni)r dr 



(C^2 _ 2^2) 

for the gaseous isopedic MSID. 

For small 7^ 0, we readily obtain 



+ 



dr2 



(86) 



(87) 



M + 



(C^2 _ 202) 



-4^2 



d 



-I- ■ 



d2 



(tu2 - 2^2 )r dr dr2 



= 
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and 



ys 



d 



+ 



_(i^2 - 2^2 )r dr ' dr2 

by removing uj ^ Q factor in equations 1861 and 1871 1. 

For the onset of marginal instability in the hmit of uj 
0, we then derive from equations 1881 and 1891 

s So (2 d dp 
^ 2Q.1\r dr dr'^ 

and 



(89) 



= 



(90) 



2ni 



2d_ ^ _<f_ 

r dr dr'^ 



$9 







(91) 



These two resuhing equations are exactly the same as 
those obtained by setting m = in equations 1571 and 1581 . 



-11 + 



!I^_2--r— ^^^-0 (92) 
(m''~2)Q.lr\ r dr dr^ ^ ' 



for the stellar SID, and 



^0 



m ^_d 



{m?-2)^lr\r " dr ' -'-^ ' ^ ° ^^^^ 

for the gaseous isopedic MSID, respectively. This situation 
of our composite MSID system parallels the case of Shu et 
al. (2000) for a single SID which is isopedically magnetized 
as well as the case of Lou & Shen (2003) for a composite 
SID system. 

For a background MSID with a coplanar magnetic field, 
the two different limiting procedures of m = and uj —> 
give rise to two different sets of resulting differential equa- 
tions (Lou 2002; Lou & Zou 2004; Shen, Liu & Lou 2005). 



4 ALIGNED AND UNALIGNED CASES 

4.1 Aligned Global Perturbation Configurations 

By definitions, aligned MHD disturbances in a disc mean 
that maxima and minima of perturbations occurring at dif- 
ferent radii line up in the azimuth. In contrast, when maxima 
and minima of perturbations in a disc shift systematically 
in the azimuthal angle from one radial location to the next, 
the resulting configuration for such disturbances are referred 
to as unaligned or spiral patterns (Kalnajs 1973; Binney & 
Tremaine 1987; Shu et al. 2000; Lou & Shen 2003; Lou & 
Zou 2004). 

4-1.1 Axisymmetric Disturbances (m — 0) 

With Lu = m = 0, the solution to equations 1411 — 1481 
takes the form of f/= = C/9 = 0, fi' = Kf/r, fi^ = Kf/r, 
J" = K^r, J9 = Kir, V" = Kllnr, and = K§lnr, 
where K^'^ for i = 1, 2, 3 are all constants which can be cho- 
sen properly such that equations H42^ . 145^ . H47|l . H48^ are 
satisfied. Nevertheless, such a solution merely represents a 
rescaling of one axisymmetric equilibrium to a neighbour- 
ing axisymmetric equilibrium. Such rescaling solutions are 
allowed but are not of particular interest here. 



4-1.2 Cases with \m\ > 1: Nonaxisymmetric Disturbances 

In discs with power-law distributions of background physi- 
cal variables, we consider aligned perturbation variables that 
carry the same power-law dependence as those of the equi- 
librium SID do. By this choice, we make use of the following 
exact potential-density pair relation, namely 



: as/r 



■ cTg/r , 



(94) 



V = -2TTGrfi' /\m\ , V = -2nGrny\m\ , (95) 

where CTs and ag are two small constants for perturbations. 
A substitution of equations 1941 and 1951 into equations 159|l 
and 1601 then leads to the following equations: 

I {Hirfi' + Gir^i") , (96) 



(i/ar^^ + Gar/i'') , (97) 



where the four coefficients H\, H2, Gi and G2 are defined by 







r 


r 


dr 


dr'^ 


rrP 




r 


r 


dr 


dr'^ 



Hi 



H2 



1 



2) 



(Dl + 1) 1 



D<>{m^^2) 

Gl EE 

G2 = - 



(1 + 5) 
e5 



{Dl + 1) 



{l + e5) 
1 



D|(m2-2)|m| {1 + 5) 

{Dl + 1) 5 
'D2(m2-2)|m| {l + eS) 



(98) 
(99) 
(100) 
(101) 



By substituting the forms of and ^® given by equation 
1941 into equations 1961 and 1971 . we readily obtain 



(1 



(102) 



(1 - JJam^)^^ = G2mV • (103) 

For nontrivial solutions of /j," and , we then derive 

(1 - Him^)(l - Ham^) = GiGam** (104) 

for the stationary dispersion relation to determined and 
thus Dg- Using definitions 1981 — 11011 . the solution condi- 
tion 11041 can be written explicitly in the form of 

+ 1 1 



Dl{rrP 



■2) 



■n\ 1 + 5 

D'g + 1 



eS 



1 + eS 



(I>?+1)(D^ + 1)771^5 



(105) 



(l + <5)(l + e<5) 
In terms of background profiles 1271 and 1281 . we obtain 

a^(l + e(5)(l + !>,') =eag(l+(5)(l + Dg) (106) 
and then 

(107) 



2^^(l + e5) a 



where the ratio (3 = a^/ag was introduced earlier. As is 
obvious, equation H105^ remains unchanged by replacing m 
with — m. For simplicity, we use m to represent |m|. In fact, 
positive and negative m correspond to trailing and leading 
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spiral patterns relative to the sense of SID rotation, respec- 
tively. After straightforward manipulations and rearrange- 
ments, equation HUSH can be cast into a quadratic equation 
in terms oi y = D^, namely 

Cay' + Ciy + Co = , (108) 

where the three coefficients C2, Ci and Co are 

/3, 



C 



e 



(|m| +2)(|m| - 1) 



2)^ + 



\m\S{e^l) 



(5 + 1) {S + iy 



Ci = 



2{\rn\ - 1) 



(m^ - 25 - 2)(1 + eS) + \m\S{e - 2 - eS) 



{1 + 5Y 

2(m^ - l)(m^ + Sm^ + |m| - 2 - 25) 



(l + '5) 

Co — — (m' — |m|) 

m^(l + 5){l + e5) + |m|(2e5^ + 5 + e5) - 2(1 -^ e5) /3 



(1-f 5)2 
2|m|(m2 - l)(|m| + \m\5 - 1) 



e 



(1+5) 

respectively. From these three coefficient expressions, the re- 
sulting y solutions remain valid for both positive and nega- 
tive signs of m. Unless otherwise stated, we shall take m to 
be positive without loss of generality. 

As seen from our MHD perturbation analysis so far, the 
Q parameter always appears together with the /3 parameter; 
it is natural and convenient to introduce a new parameter 
X = /3/6 to simplify notations. Physically, Qa^ represents a 
sort of magnetosonic speed squared given the isopedic mag- 
netic field geometry; the ratio x = 0.3 /{^oig) thus stands for 
the square of the ratio between the stellar velocity dispersion 
and the magnetosonic speed. 

We now discuss properties of quadratic equation IIU8II 
with its determinant A given explicitly by 

A = Ci-4CoC2-^^^-^ 

x{[(m + 2)(m-l)(xe5-l) 

+ -2){x-S)Y +4m^xS} ■ (109) 

Except for the case of m = 1, the important fact that A > 
for m > 1 means that quadratic equation IIU8I I always have 
two real solutions of y = D^. The physical meaning is that 
for aligned nonaxisymmetric stationary configurations to ex- 
ist in a composite MSID system with an isopedic magnetic 
field, condition 11U5L or equivalently, condition 111)811 must 
be satisfied for appropriate values of given a set of spec- 
ified parameters m, 5, e and x = The two real y = 
solutions may not be necessarily physical unless the nonneg- 
ative requirements of both > and Dg > are met. 

We now elaborate several subtle points below. 

For the aligned case with m = 1, it turns out that C2 = 
Ci = Co = by definitions. Equation l|108|l can therefore be 
satisfied for arbitrary values of D^. This situation is quite 
similar to the m = 1 cases of a single isopedically magnetized 



SID studied by Shu et al. (2000) and of a composite system 
of two coupled SIDs analyzed by Lou & Shen (2003). 

When m > 2, the expressions of coefficient parameters 
in the equation can be rearranged into the following forms 
by multiplying through a factor of (1 + 5)^, namely 

C2 = x{m + 2){(m' - 2)(65 + 1)(5 + 1) 

+m[(5e + l)-(5 + l)]} , 
Ci =2x{[m'-2(5 + l)](l + e5) 

+m[2(5e + 1) - (5 + 1) - (e5 + 1)(5 + 1)]} 
-2(m+l)(5 + l)[(5 + l)(m' -2) + m] , 
Co = ~mx{m\l + 5)(1 + e5) - 2(1 + e5) 

+m[2(e5 + 1)(5 + 1) - (5 + 1) - (e5 + 1)]} 
+2m(5 + l)(m + l)[m(5 + 1) - 1] . 

Since 1 + 5 >0, l + e5 >0 and x > 0, we further simplify 
these three coefficient expressions to obtain equivalently 



C2 



Ci 



Co 



(m + 2) 



- 2 



5 + 1 
2(m + 1) 



(5e + l) (5 + 1) 
2 



1 



i^S + l)x 
-m m' + m ( 2 
2m(m+ 1) 



5 + 1 
[(5 + l)(m" - 
1 



5 + 1 
m(5 + 1) - 



5e + 1 

2)- 



e5 + 1 



- 1 



iS + l) 



1] 



(110) 



(111) 



(112) 



{^S + l)x 

For the convenience of further analysis, we introduce two 
new dimensionless parameters B and X defined by 

(1 + ^) 



X 



- (l + eS) 
m{B-l) 



(1 + 5) 

It then follows that C2, Ci and Co can be expressed as 
C2 = (m + 2)X , 

C.^2X-2(^ + l)[^+^^ 

1 



113) 
114) 



- 2 + 



Co = -m\ X -2{m + 1) 



B 

— [m — 
X\ 1 + 5 



1 + 5 
5 



(1 + 5) 

Because there are several parameters in our theoreti- 
cal model analysis and our main motivation is for galactic 
applications, it would be efficient and more sensible to have 
rough ranges to bracket these parameters in order to explore 
different regimes for various possible solutions numerically. 

For this purpose, we estimate ranges of the parame- 
ters A, (3, e and for disc galaxies. A typical magnetic 
field in a spiral galaxy is about 1 ~ 10/iG and higher 
values of ~ 30 — 40/iG may be reached in central cir- 
cumnuclear regions by equipartition arguments (e.g., Lou 
et al. 2001). A gas surface mass density is about ~ 
2 X 10"'' - 2 X 10"^ g cm . The value of mass-to-flux ratio 
A is about 20 ~ 2000 g cm"'G"\ Therefore, the A param- 
eter falls in the range of 0.03 ^ 3. Here, we estimate the 
maximum range for every parameter in the galactic context 
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while keeping in mind the constraint of 1 + e(5 > for the 
background magneto-rotational equilibrium. The variation 
range of e may be as large as —1110 ~ 0.9. Taking 5 as 0.01, 
we choose the e parameter to fall within the rough range of 
— 100 ~ 0.9. In a typical late-type spiral galaxy, the stellar 
velocity dispersion Os is usually several times higher than 
the sound speed ag of the gaseous disc, such that /3 > 1 in 
our model analysis. We have aheady shown the O range to 
be 1 < < 2. As O and /3 always appear together, we use a 
single combined parameter x = and further take x > 1 
in reference to/3>l. Asl<0<2, the effect of Q varia- 
tion on X is somewhat limited. For example, one may take 
ag = 7 km s~^ and as — 30 km s~^ for a typical late-type 
spiral galaxy. It is then justifiable to take x = ^^/(Qig) > 1- 
The introduction of B parameter by definition H118|l is es- 
pecially useful to estimate effects of an isopedic magnetic 
field. It is clear that inequality B > 1 holds valid given the 
condition e < 1 and the requirement l + e5 > for the back- 
ground magneto-rotational equilibrium. With a fixed value 
of 5, the increase of the isopedic magnetic field flux will lead 
to an increase of B. It turns out that Bis a, fairly good indica- 
tor for the isopedic magnetic field, although the relationship 
between the two is not as simple as a linear one. 

Formally, the two real y = solutions of quadratic 
equation l|108|l are given by 



J/1,2 



2C2 



(115) 



We now introduce an important parameter Ra (see Ap- 
pendix A for more details) defined by 

m[m- 1/(1 + 5)] 



Ra 



m2-2 + m/(l + 5) 



> 



(116) 



Based on the mathematical conclusion reached in Appendix 
A, for the yi solution with the plus sign in II115II . we have 
yi > -Ra > when C2 > and yi < -1 when C2 < 0. 
For the 1/2 solution with the minus sign in lllSII . we have 
-1< 2/2 < Ra- 
in the previous analyses (Lou & Shen 2003; Shen & 
Lou 2003, 2004a, b; Lou & Zou 2004), the two solutions yi 
and y2 are sometimes referred to as supersonic and subsonic 
rotation solutions respectively, because as rotational Mach 
numbers, yi is almost always greater than 1 and j/2 remains 
always less than 1. In comparison, the situation becomes 
somewhat different in the presence of an isopedic magnetic 
field. In particular, yi may drop below 1 for some specific 
values of m, 5, B and x\ this could also happen when there 
is no magnetic field. Furthermore, with the boosting of the 
magnetic field, y2 may become greater than 1 even under the 
condition x ^ 1- Note that 7/2 remains always less than 1 in 
the absence of magnetic field. Thus for proper terminologies, 
we use the F-wave solution for yi and the S-wave solution for 
j/2, because when they are both physical, yi and 1/2 remain 
always faster (F) and slower (S), respectively. 

When B = 1 for the absence of magnetic field, we have 
= 1 and hence x ~ P- This situation is exactly the same 
as the problem analyzed by Lou & Shen (2003) which will 
be frequently referred to in our following discussions. 



For the F-wave solution j/i, there exists an upper limit 
for the magnetic field intensity. When the magnetic flux be- 
comes so strong such that B has a sufficiently high value, C2 
will become negative and thus yi becomes unphysical for 
being negative. The physical constraint on B for the F-wave 
solution yi is therefore 

(m^ -2)(l + 5) 



1< S < 1 



(117) 



where the lower bound on the left-hand side (LHS) is from 
the definition of B while the upper bound on the right-hand 
side (RHS) is derived by requiring C2 to be positive. We note 
that with the increase of both m and S, the allowed range 
of B for the existence of an F-wave solution yi becomes 
enlarged. Moreover for yi > 0, yi decreases with increasing 
X- 

For the S-wave solution 2/2, the somewhat loose con- 
straint is that —1 < 2/2 < Ra- For a physical solution, we 
should require 2/2 > 0. Since j/i > for C2 > and 2/1 < for 
C2 < 0, we infer from the solution property j/12/2 ~ C0/C2 
that J/2 > for Co > and 2/2 < for Co < 0, respectively. 
For a physical S-wave solution 2/2 > 0, the condition Co > 
is simply equivalent to the following inequality 
'2{m+l) 



B 



+ 



- (m + 2) > 



(118) 



X m(l + 5)-l_ 

When B = 1 for the absence of magnetic field, by setting 
the LHS of the above inequality III8II equal to zero, we 
determine a critical value Xc for x which is the same as /3c 
parameter of Lou & Shen (2003); as /3 or x is increased to 
become greater than /3c or Xc the S-wave solution t/2 will 
change from positive to negative values. Furthermore, the 
increase of B leads to an increase of Xc- For B > 1, the 
critical value Xc is given explicitly by 

_ 2(m + l)[m(l + 5)-l]g 2(m + l)g 
~ {m + 2)[m(l + S) - 1] - mB (m + 2) ' ^ ' 
One interesting point to emphasize is that when B becomes 
sufBciently large, there is no limit on x and the S-solution 
2/2 remains always positive. More specifically, when the fol- 
lowing condition 112011 is fulfilled 



B>1 



{rn^ -2)(l + 5) _^ 2(m + l)5 



(120) 



m m 
inequality 111811 will be always satisfied and thus the S- 
solution 2/2 remains always positive. 

Without magnetic field, there always exists at least one 
positive solution (see Lou & Shen 2003). In contrast, when 
an isopedic magnetic field is anchored in a gaseous SID, there 
does exist a combination of parameters such that both 2/1 
and 2/2 solutions may become negative. Under these circum- 
stances, there would be no stationary MHD density wave 
patterns supported by SID rotation. For this to happen, the 
range of B parameter is given below 



^ I jm^ - 2)(1 + 5) ^ g ^ ^ I (m" - 2)(1 + S) ^ 2(m + 1)5 

m mm 

(121) 

[see the RHS of inequality l|117|l and inequality H120^ ]: the 
LHS of inequality 112H leads to 2/1 < 0, while the RHS of 
inequality 112111 makes 2/2 < possible and the x parameter 
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Figure 1. Two sets of solution curves for y = versus x varia- 
tion for the aligned case with |m| = 2 and B = I for the absence 
of magnetic field. We use B in the figure to correspond to B in 
the main text. The vertical line marks the constraint of x > 1- 
We take 5 = 0.2 (solid lines), 1 (dotted lines) and 5 (dash-dotted 
lines), respectively. The allowed range for x [li+oo)i a^nd 
plot these curves within the x interval (0,20]. This figure is the 
same as figure 1 of Lou & Shen (2003) and serves as our reference 
to understand the role and effects of an isopedic magnetic field. 



should be adjusted to make Co < such that y2 < 0. We 
take parameters m = 2,S = l,B — 4 and x = 20 as a specific 
example, and the two negative solutions are yi = —2.8651 
and j/2 = —0.3490, respectively. 

For the phase relationship between the two surface mass 
perturbations /i" and fi^ , we note that equation ll()2|l can 
be written in the form of 



= -1 - 



= -1 - 



[Dlim^ - 2) - m\l + 5) 



|m|(Di + l) 
2)(l + (5) 2(m2-l)(l 



-5) 



\m\ |mj(i32 + i) 

2)(l + J) + |m|] {D^-Rg) 



(122) 



|m| Pi + 1) ■ 

According to expression l|122|l of fi^ /jj," , it is fairly clear that 
for the F-wave solution (i.e., yi > Ra with the unphysical 
case of j/i < being ignored), the ratio fi^ /fi" remains always 
negative, while for the S-wave solution j/2, the ratio A'^/m^ 
remains always positive. Physically, when the two density 
perturbations are out of phase, the gravity effect is weaker 
and the MHD density wave speed is faster, while when the 
two density perturbations are in phase, the gravity effect is 
stronger and the MHD density wave speed is slower (Lou & 
Fan 1998b). 

In addition to the requirement of > 0, we must also 
make sure of Dg > such that the relevant MHD density 
wave mode is physically plausible. We now write equation 
110711 in a simple form of 



(123) 



Without magnetic field with B = 1 and under the assump- 
tion of X = /3 > 1, it is easy to show that Dg remains always 
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Figure 2. Two sets of solution curves oi y = u% versus x varia- 
tion for the aligned case with |m| = 2 and f? = 4 for the presence 
of an isopedic magnetic field. We take <5 = 0.2 (solid lines), 1 (dot- 
ted lines) and 5 (dash-dotted lines), respectively. Again, we use 
a in the figure to correspond to B in the main text. The vertical 
line marks the constraint of x > 1- Note that the F-wave solution 
y\ with 5 = 1 is not shown as it drops far below zero. In order 
to give an intuitive comparison between Fig. 1 and Fig. 2 here, 
we use the same scale for the reference frame. It should be noted 
that only a range of x is permitted when the S-wave solution j/2 
is positive due to the constraint of Xa ■ 



positive according to equation H123^ (Lou & Shen 2003; Lou 
& Zou 2004). However, this may not be true in the pres- 
ence of an isopedic magnetic field with B > 1 and Dg may 
become negative even for x > 1- In other words, one more 
constraint of Dg > should be checked carefully. On the 
basis of Appendix A, we here define a new critical value Xa 



Xa 



B 



2(m + 1) 



2 - 



mS 



B 



(124) 



(m-l)(l + 5) 

and only when X > Xa can inequality 7?^ > be guar- 
anteed. In the absence of magnetic field (i.e. B = 1), we 
readily see that Xs < 1. So a more general condition for a 
physical solution would be x > niax{l,X3}. This limit of 
Xa emerges as a constraint in the presence of an isopedic 
magnetic field (B > 1). With the increase of the magnetic 
flux (i.e. B grows), the value of Xs increases accordingly. By 
comparing the expressions of Xc and Xs, it should be noted 
that wherever there exists an upper limit Xc this Xc remains 
always larger than Xa- 

We now review several properties of the two solutions 
of y = D^ as a summary. In the absence of magnetic field 
with B = 1 (see Figure 1), our analysis here is the same 
as those of Lou & Shen (2003). The F-wave solution yi re- 
mains always positive and the S-wave solution j/2 changes 
from positive to negative as x becomes larger than Xc In 
the presence of an isopedic magnetic field with S > 1, we 
define a Xs instead of Xc For a physical configuration solu- 
tion, we then require x > max{l, Xa}- With the increase of 
the magnetic flux (i.e., an increase of B), the two solutions 
2/1 and y^ become larger and the value of Xc increases ac- 
cordingly. However, as B exceeds a certain given value [see 
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expression Ill7ll 1. the F-wave solution j/i becomes negative 
and thus unphysical (see Fig. 2). In addition to the require- 
ment 1 + e5 > on the magnetic field for the background 
magneto-rotational equilibrium, a further condition 1)117^ is 
needed for the existence of an F-wave solution j/i. As Xs 
becomes larger than 1, the requirement of x > rna^xfljXs} 
means that only if Xs < X < Xc, can the S-wave solution 
j/2 be physically valid. Note that these constraints on x do 
not apply for the F-wave solution yi. As the magnetic field 
strength is increased further [see equation (I12UII 1. the up- 
per limit Xc for x will disappear and only the lower limit 
Xs exists (see Fig. 3). It is generally true that when the 
two solutions j/i and 3/2 are both positive and physical, they 
both increase with the increase of B and decrease with the 
increase of x- 

Physically, for the F-wave solution with fi^ / < 0, the 
surface mass density perturbations of stellar and gaseous 
SIDs are completely out of phase to reduce the effect of 
gravity. Meanwhile, a faster azimuthal density wave speed 
corresponds to a faster SID rotation (larger D^) in order to 
maintain a stationary MHD configuration. For the S-wave 
solution with ^j.^ / jJ." > 0, the surface mass perturbations in 
gaseous and stellar discs are in phase with each other. As 
the effect of self-gravity is enhanced, the azimuthal density 
wave speed becomes slower corresponding to a slower SID 
rotation (smaller D^) in order to sustain a stationary MHD 
perturbation configuration. 

The isopedic magnetic field plays an interesting role. 
When the magnetic field grows stronger as B becomes larger, 
the F-wave solution yi and the S-wave solution j/2 both 
become larger. According to expression 11221 for the ratio 
fj,^ / , it is easy to demonstrate that fJ,^ / i-i" will decrease 
with the increase of for both j/i and 3/2 solutions. For 
the F-wave solution, the decrease of ^j.^ / fi" means an in- 
crease of the absolute value of jj,^ / fi" , while for the S-wave 
solution, this ratio magnitude becomes smaller. A physical 
interpretation is that for a given composite SID system, an 
isopedic magnetic field will amplify the contrast of surface 
mass density perturbations in the two SIDs for an F-wave 
mode and will reduce such contrast of surface mass density 
perturbations in the two SIDs for a S-wave mode. 

As noted earlier (Lou 2002; Lou & Fan 2002; Lou & 
Shen 2003; Lou & Zou 2004), we here emphasize again the 
physical perspective that stationary aligned MHD perturba- 
tion configurations should be regarded as purely azimuthal 
propagation of MHD density waves counterbalanced by the 
advection of (M)SID differential rotation. Together with 
equation 1107II . we may write equation 11051 in the following 
form to derive the two y solutions. 

("^'-2+^)(0?-i?.) 

\m + 2 I + 5J \ 1 + 5/. 

^{Dl + l){Dl + l)^^ , (125) 

where C2 and Ra are defined by expressions IllOt and III6II , 



respectively. For the F-wave solution t/i, we have > Ra 
when > which is guaranteed by C2 > 0. Moreover from 
equations 112511 and 11231 1. we can show that for a positive 
F-wave solution j/i, the corresponding Dg is also positive 
such that the overall solution is physically plausible. 

We now write solution condition 11051 for stationary 
MHD configurations of a composite SID system in a physi- 
cally more suggestive form of 

[n'iim'^ - 2) - m^a^/r^ + 2-!vGT,o\m\/r] 

X [J7g(m^ - 2) - m'^eal/r'^ + 27reGEg |m|/r-] 

= 47r'G^E^EgmVr' . (126) 

The RHS of equation 112611 represents the mutual gravita- 
tional coupling between the 'fiuid' stellar SID and the 'mag- 
netofiuid' of a gaseous isopedic MSID. Without this gravi- 
tational coupling, the two factors on the LHS would emerge 
as two separate solution conditions for stationary aligned 
perturbation configurations with |m| > 2 of stellar SID and 
of gaseous isopedic MSID, respectively. For the stellar SID 
alone, the condition would be 

m^Q,^ = hii + m^a^/r^ — 27rGEo|m|/r- (127) 

and for the gaseous MSID alone with an isopedic magnetic 
field, the condition would be 

m^nl = Kg + m^Qal/r^ - 27reGEg|mI/r . (128) 

According to the well-known dispersion relation of density 
waves first derived under the tight-winding or WKBJ ap- 
proximation [Lin & Shu 1964, 1966; Lin 1987; or equation 
(39) of Shu et al. 2000], dispersion relations lfT?7^ and l(T^ 
can be readily recovered by replacing the radial wavenumber 
fc| with the azimuthal wavenumber \m\/r and setting u — 
in an inertial frame of reference as noted by Lou (2002) in the 
study of stationary MHD perturbation configurations of a 
single MSID with a coplanar magnetic field. Along this line 
of reasoning, we know that dispersion relations 112711 and 
112811 separately describe azimuthal propagations of hydro- 
dynamic density waves in a stellar SID and of MHD density 
waves in a gaseous MSID with an isopedic magnetic field. 

4- 1-3 Secular Barlike Instabilities 

We have constructed analytically stationary configurations 
for aligned nonaxisymmetric MHD perturbations in a com- 
posite MSID system. Whether these perturbation configu- 
rations are stable or just represent transition states from 
axisymmetric equilibria to non-axisymmetric configurations 
still remains an open question (see discussions of Shu et al. 
2000 and extensive references therein). 

It was tentatively suggested (Shu et al. 2000; Galli et al. 
2001) that in the case of a single (M)SID, these stationary 
solutions signal onsets of bifurcations from an axisymmetric 
(M)SID to non-axisymmetric configurations in (M)SIDs. 

Following the earlier procedures (Shu et al. 2000; Lou 
2002; Lou & Shen 2003), we explore parameter regimes of 
these stationarity conditions in reference to the disc stabil- 
ity criterion first hypothesized by Ostriker & Peebles (1973) 
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Figure 3. Two sets of solution curves for versus B varia- 
tion. Here we use B in tlie figure to correspond to B in the main 
text. The vertical line represents the diverging point of yi for the 
aligned case with |m| = 2, (5 = 1 and x = li 5 and 10, respectively. 



for the onsets of bar-type instabilities. Judging similarities 
and differences between a single (M)SID and a composite 
(M)SID system, we examine stability properties for station- 
ary configurations of aligned and unaligned perturbations 
with an isopedic magnetic field in reference to the analysis 
of Shuet al. (2000). 

The criteria for (secular and dynamic) barlike insta- 
bilities were sometimes expressed in terms of the ratio of 
the rotational kinetic energy T to the absolute value of the 
gravitational potential energy W (Ostriker & Peebles 1973; 
Binney & Tremaine 1987; Shu et al. 2000; Lou & Shen 2003; 
Lou & Zou 2004). To obtain the scalar virial theorem in the 
present case of a composite MSID system, we start from the 
radial force-balance equations of the equilibrium state, 



E=J7?r=^(a?E;^) + SS^( 



Sgf^'-=i(e«^Sg)-fSgJ-(<^^+e0«) 



(129) 



(130) 



where the axisymmetric dark matter halo potential is not 
included yet. An addition of equations I129II and 113UII gives 

(Egn^ + E2J7?)r = ^(ea^Eg + a^Eg) 

+^±{cb' + ecbn + ^o^(<P' + ^n ■ (131) 

Multiplying equation l|131|l by a ring area element 2Tvr^dr 
and integrating from to a finite radius L, we obtain 

2(T + W) + W = 27rL^[eagEg(L) +a?Ef)(L)] , (132) 



where 



T: 



]^T.l{rQ.gf2-Krdr + 



(ea^Eg -I- ajEo)27rrdr , 



^2-Krdr , (133) 
(134) 



r-E 



dr 



2Tirdr 



rEi 



d(0" + 
dr 



■2'Krdr 



(135) 



Here T is the rotational kinetic energy in the composite SID 
system, U is the equivalent 'thermal energy' contained in 
the composite MSID system including the effect of magnetic 
pressure, and W is the equivalent gravitational- work integral 
including the effect of magnetic tension. 

Using S = Eg/Eg, fls = asDs/r, = oiagDg/r to- 
gether with equations I27II and 12811 . we cast these integrals 
in the following forms of 

al(l + D?) 



T 



U - 



2G{l + S) 



Di 



^(i±^(l + D^) 

{1 + 5) + 



95 



Q5_ 

P 



1 L 



G(l + (5)2 



-(l + 2(5 + er)L 



(136) 
(137) 
(138) 



For a composite MSID system of an infinite radial ex- 
tent, all three integrals 11361 — 1113811 above diverge as L ^ 
+CXD but their mutual ratios remain finite. For example, the 
ratio of the kinetic energy of disc rotation to the absolute 
value of the gravitational potential energy is 



T 

m 



(1 + Dl){l + 25 + 6^2) _ (1 + 5)[i + [Q5/I3)] 



2(l-f Di)(l + 25 + e52) 
1 + {5/x) 



(139) 



2 2(l + D2)(l + 5/B) ■ 
This ratio can also be arranged into the following form 
T _1 al-^iDl + QalnDl 

|W| 2[aiEg(l + Z?i) + eaiEg(l + D|)] ^ 
which is explicitly symmetrized with respect to physical pa- 
rameters of the two SIDs. Here, e does not appear explicitly 
in the final expression. 

Note that the ratio T/|W| falls between and 0.5 as 
usual (Binney & Tremaine 1987; Lou & Shen 2003; Lou & 
Zou 2004) and increases with increasing D^. For stationary 
configurations of aligned MHD perturbations in a composite 
system of two coupled SIDs, the two possible values of 
give rise to two different values of T/|W| ratio; the larger 
and smaller values of correspond to larger and smaller 
values of the T/|W| ratio, respectively. 

Furthermore, as the magnetic flux increases (i.e., an in- 
crease of B), the T/|W| ratio decreases. This simply means 
that the effect of an isopedic magnetic field tends to make 
the composite SID system more stable (Lou & Shen 2003; 
Shen & Lou 2003, 2004a, b; Lou k Ion 2004, 2005). 

On the basis of an extensive numerical exploration for 
m = 2 to m = 10 and so forth, one can show the exis- 
tence of a dividing line in the ratio T/|W| for global sta- 
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Figure 4. The T/|>V| ratio versus x variation for the aligned 
cases with m = 2, B = 1 and 5 = 0.2 (solid lines), 1 (dotted lines) 
and 5 (dash-dotted lines), respectively. Here, we use B, W and 
T in the figure to correspond to B, W and T in the main text, 
respectively. The horizontal line raj (Am + 4) = 0.1667 divides 
the ratios T/\W\ for the S-wave solutions y2 and for the F-wave 
solutions yi. Through numerical explorations for m = 2 up to 
m = 10 and so forth, the ratio T/|W| of the S-wave solution j/2 
reaches a maximum value of m/(4m + 4) at x = 



tionary configurations of aligned MHD perturbations. More 
specifically, for F-wave solutions yi, the ratio T/|W| ranges 
from m/(4m + 4) to 1/2, while for S-wave solutions y2, the 
ratio T/|>V| ranges from to m/(4m -I- 4) when m > 2. 
Based on A'^— body numerical simulation experiments in- 
volving only ~ 300 particles for a stability analysis of a 
rotating disc under self-gravity, Ostriker & Peebles (1973) 
suggested an empirical criterion that T/|W| < 0.14 ± 0.02 
is necessary but not sufficient against bar-type instabilities 
(see also earlier numerical simulation results of Miller et al. 
1970; Hohl 1971). When T/\W\ > 0.14 ± 0.02, a disc sys- 
tem would evolve rapidly into bar-type configurations (e.g. 
Binney & Tremaine 1987). In the current problem of a com- 
posite MSID system, the lowest value of T/|VV| ratio for 
F-wave solutions yi is ~ 0.1667, while for S-wave solutions 
3/2, the T /\W\ ratio can be lower than 0.14 ± 0.02. 

In the above analysis, we have not included the effect 
of a massive dark matter halo. One obvious consequence 
of a massive dark matter halo is to increase \W\ and thus 
decrease T/|W| ratio, and for both classes of global pertur- 
bation configurations, this implies the tendency towards a 
stability. At any rate, a comparison of the above two cri- 
teria suggests that given a background dark matter halo, 
F-wave configurations tend to be more stable than S-wave 
configurations. 

By numerical experiments for m = 2 up to m = 10 and 
so forth, we found that when B — Xy the value of T/\W\ 
ratio for S-wave solutions 7/2 is always m/(4m + 4). 



4.2 Unaligned Logarithmic Spiral Disturbances 

For global stationary configurations of unaligned or loga- 
rithmic spiral MHD disturbances, we take on the following 
set of exact density-potential pair relations that satisfy the 
Poisson integrals, namely 



-3/2 



exp(ia In r) 



(141) 



V = V r 



"^''^exp(ialnr) , (142) 

where a' and v' are small constants for perturbations, q is a 
parameter to characterize the radial variation; the two small 
coefficients and a' are algebraically related by 

v' = -27rGA/'™(a)cr' , (143) 

with J^rn(ci) = K{a, m) being the Kalnajs function (Kalnajs 
1971) and superscript I = s,g. For logarithmic spiral pertur- 
bations in a composite MSID system, the above description 
is a sensible extension of earlier model studies. The radial 
scaling parameter a (closely related to the radial wavenum- 
ber) is naturally taken to be the same for perturbations in 
both stellar SID and gaseous isopedic MSID. 

In the following analysis and computations, we use two 
useful formulae of J\fm{ct). One is the recursion relation in 
m of Afm{a) for a fixed a value (Kalnajs 1971) 



AA™+i(Q)AA„(a) = [(m + 1/2)" + 



(144) 



and the other is the asymptotic expression of J\fm(a) 

AA™(a)Ri(mVa' + l/4)-^/2 (145) 

in the regime of -f q2 > 1 (e.g., Shu et al. 2000). When 
the requirement for accuracy is not so stringent in some 
quantitative analyses, expression 114511 may even be used to 
compute logarithmic spiral solutions with \m\ — 1. 

For |m| > 0, we now proceed to solve equations 1591 
and I6OI I that can be cast into more compact forms of 

'(i/ir/i= + GiV) , (146) 



m 
r 



dr 



dr'^ 



where the four coefficients H\, H2, Gi and G2 are 



Hi 



H2 



Di(m2 -2) 
1 



(75? + l)AA„(a) 



1 - 



(1 + 5) 
(D^ + l)e5jVm{a) 



Gi 



G2 



(l + eS) 
D1{m^ - 2) (I + 5) 
Dj{m^ - 2) (1 + eS) 



(147) 

(148) 
(149) 
(150) 
(151) 



Substituting jj." and /i^ in the form of expression 114111 into 
equations 114611 and 114711 . we readily obtain 

[1 - Hi{m'^ +a'^ + l/4)]At" = Gi(m^ + + l/4)^t^ , (152) 
[1 - H2im + l/4)]/i» = G2(m^ + a + l/4)^t" . (153) 

Multiplying both sides of the above two equations and re- 
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moving fi^i-t^ 7^ 0, we derive the solution condition as 

[1 - Hi{rn +a+ 1/4)][1 - H^im^ + a+ 1/4)] 

= GiG2(m2 + a' + l/4)2 . (154) 

Using definitions H148^ — Ijl51|l . condition Ijl54|l then becomes 

1 



I>2(m2 - 2) 



I>2(m2 - 2) 



_ (Og + 1)A/-^(q) 
(1 + 5) 



1 - 



(l + eS) 



2 , 2,1 

4 



2 , 2,1 

m + Q H — 
4 



2 , 2,1 

m + a + - 
4 



DiD2(m2-2)2 (l + 5)(l + e5) 

This equation may be rearranged into the form of 
r / n2 



(m'' + 1/4) I 



(D^ + l)e<5AA™(a) 



= — {i+m + es) — +^/^) 

Using expression HUTU for Dg, we readily obtain a quadratic 
equation in terms of y = D^, 

Cay' + Ciy + Co = , (155) 

where the three coefficients C2, Ci and Co are 



,2 r,^ , J^m{a) ( 2 , 2,1 



(l + e<5) , 2 , e5Mm(a) ( 2 , 2,1 

'?AC.(q) ^ 2^ 2^ 1 

(1 + 5)2V 4 



Ci =X 



(1 + 5) (1 + 5) Iv 4 

2Mn.{o)( 2^ 2^ l^ 2 , 9 



(1 + 5) 



25Ae(Q) 2^ 2^ 1 

iTT^Fr +4 



/ 2 o\ , -'^m(") 2 , 2,1 



2m' + a' — 



Co = X ( 'Ti' + a' + i 



A/"m(a) 
(l + <5) 



(1 + 5)^ ^ 

e<5Mn(«) ( 2 , 2,1 
/ 2^ 2^1' 

(Twr +4 



2m +a m +a H — 

4 A 4 



A^m(Q) _ ' 
(1+5) 



In the following, we shall remove some nonzero common fac- 
tors of C2, Ci and Co to simplify these expressions. For this 
purpose, we introduce a new parameter 

(m^ - 2) 



M 



{m^+a^ + 1/4) 



(156) 



Together with parameter B defined by equation 111311 as in 
the aligned case, we obtain 



Ci 

'(1 + 5) 
Co 



C2 = [M+Mm{a)\ 
1 



B 



\M+Nn,{a)] 



M 



[M + 2Nm{a)-l] \ - 



1 + 5 
2Mm{a 
(1 + 5) 
(A4 + 1) 



A/'m(a) 



M + 



(1 + 5) 



(1 + 5) 



■\M+Um{a)\ 
7^[AA™(«)-1] 



(1 + 5) 
{M + 1) 



Nm{a) 
(1 + 5) 



- 1 



Consequently, we have the determinant A of quadratic equa- 
tion p55|l above as 



Nm{a) 
"(1 + 5) 



4C0C2 = {M + lf^[^[M+Nm{a)\ 



M + 



(1 + 5) 



2^ 45Ae(a) 



(l+5)2x 



> 



Formally, there are two different real solutions j/i and 2/2 of 
quadratic equation H155II in the forms of 

-Ci±(C?-4CoC2)i/' 



yi,2 



2C2 



(157) 



for A > 0. 

Based on our earlier discussion for the axisymmetric sit- 
uation of m = in equations 188^ and H89|l for both aligned 
and unaligned cases, the above results remain valid for the 
axisymmetric case of m = in the unaligned case (i.e. per- 
turbations with radial oscillations). Note that A remains 
always positive for m > 1. The rare case of A = can only 
happen when m = Q and + 1 = 0. In general, there are 
therefore two different real solutions y\ and 1/2 and these two 
solutions may become the same only if m = and A^ + 1 = 0. 

We now proceed to discuss the three situations m > 2, 
m = 1 and m = in three separate subsections below. 



^.2.1 Logarithmic Spiral Cases with m > 2 

When m > 2, we take the following approximate form of 
7Vm(a), namely 

•^'"(") - („^2 + q2 + 1/4)1/2 • 

As in the aligned case, we introduce a similar parameter 



Rs = 



[1-A/'^(q)/(1 + <5)] 



(158) 



(see Appendix B for more details). 

On the basis of the mathematical analysis in Appendix 
B, for the j/i solution, we have yi > Rs > Q when C2 > 
while j/i < — 1 when C2 < 0, while for the 1/2 solution, we 
have — 1 < j/2 < -Rs independent of the sign of C2. Here, the 
role of this Rs parameter for the unaligned case is similar to 
the role of Ra parameter for the aligned case. Note that when 
m > 2, parameter Rs remains always positive. As in the 
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aligned case, we still refer to y\ as the F-wave solution and 
2/2 as the S-wave solution for the unaligned case, respectively. 

To analyze the phase relationship between the two sur- 
face mass perturbations ^j," and in the coupled stellar SID 
and isopedic MSID, we recast equation 11152^ in the form of 



[l^Jj-i(m^+a" + l/4)] 
Gi(m2 + q2 + 1/4) 

M{l + 5) {M + l){l + 5) 
Um{a) Mm{a){D-i + l) 
\Nn,{a) + M{1 + 5)] (D'i-Rs) 



(D'i + 1) 



(159) 
(160) 



From expression 11601 of ratio jj,^ / fi" , it is clear that for 
the F-wave solution yi (i.e., yi > Rs with the unphysical 
situation of yi < being ignored), ratio n^/fi" is always 
negative, while for the S-wave solution y2, ratio fi^/fi" is 
always positive. Note that these results are qualitatively 
the same as those in the aligned case with m > 2. The 
same as the aligned case, for the F-wave solution j/i, ratio 
IJ.^ /iJ." < means that surface mass density perturbations in 
stellar SID and gaseous MSID are out of phase while for the 
S-wave solution y2, ratio fi^/fi" > means that the surface 
mass density disturbance in the gaseous MSID is in phase 
with the surface mass density disturbance in the stellar SID. 
It becomes transparent that in-phase and out-of-phase sur- 
face mass density perturbations are intrinsic characters for 
the two solutions yi and y2 when m > 2. Physically, both 
aligned and unaligned perturbations represent MHD density 
wave propagations in a composite MSID system. As MHD 
density wave patterns with out-of-phase surface mass den- 
sity perturbations travel faster, a higher rotational Mach 
number Ds is required in order to strike a global stationary 
configuration. In contrast, MHD density wave patterns with 
in-phase surface mass density perturbations travel slower, a 
lower rotational Mach number Ds is needed in order to strike 
a global stationary configuration. This phase relationship for 
the two surface mass density perturbations is retained even 
in the presence of an isopedic magnetic field (Lou & Fan 
1998; Lou & Shen 2003; Lou & Zou 2004). 

The solution y stands for which is the square of 
the rotational Mach number Ds of the stellar SID. The two 
solutions yi and j/2 correspond to the square of two different 
rotational Mach number. One obvious physical requirement 
is that both yi and j/2 must be positive. Meanwhile, we also 
need to require Dg to be positive in order to establish a 
plausible composite MSID system. According to equation 
II123II . we have the correspondence between the condition 
Dg > and the condition 



Ds = y> — 
X 



(161) 



We discuss in the following these two constraints for F-wave 
solutions yi and S-wave solutions j/2, separately. The two 
constraints must be satisfied simultaneously to limit the pa- 
rameter space, otherwise there are no physical y solutions. 
For F-wave solutions j/i to exist, we have the following 



mathematical conclusions from Appendix B, namely 



yi>0 



and 



yi > 



X 



are guaranteed by C2 > 0. As a result, we can set a limit 
on the strength of an isopedic magnetic field by requiring 
C2 > 0. For the B parameter, we thus have inequalities 

M{l + S) 



1 < B < 1 -I- 



(162) 



A/'m(Q) 

for the existence of a physical F-wave solution j/i . As in the 
aligned case, this constraint means that too strong an iso- 
pedic magnetic field would be impossible for sustaining a 
stationary F-wave solution yi. By increasing both m and 
S and decreasing a, the allowed range of B for sustaining 
F-wave solution j/i becomes enlarged. Once the F-wave so- 
lution yi exists, we always have a decreasing yi as both x 
and B increase (see Appendix B for details). Through our 
analysis, we also find that for B > 1 and q 00, the C2 
coefficient will become negative and thus lead to a negative 
2/1 . Therefore, there exists a critical value Oc of a such that 
2/1 > exists only when a < ac; this critical Oc is 



2f{l + Sf 



2 I 
m + 



1 1/2 



(B-l)2 V"^ ' 4, 

For the S-wave solution 2/2 > to exist, we establish the 
following two correspondences between 



y2 > 
and between 



2/2 > 



B 

X 



and 



C2(^ 



Co > , 



and 



+ Ci(^ 

X 



1 + Co > 



respectively. The requirement of Co > for a positive 2/2 
leads to a constraint on x parameter, that is, x < Xc where 
the critical value Xc for x is explicitly defined by 

BjM + 1)[1 + S - M„,{a)] 

[l + 5-J\fm (a)] [M + (a)] - SAA™ (a) [1 - (a)] 
B{M + 1) 



Xc 



We further note that when B parameter becomes sufficiently 
large to satisfy the following inequality 

B>l + :^lp^ + #±^, (164) 

the X parameter will be free from the constraint of x < Xc 
because such Xc no longer exists. 
From the inequality 



(163) 



C2 



we derive another constraint on x, that is, X > Xs where Xs 
is explicitly defined by 

B 



Xs 



{M + l) 



M+Nm[a) 



(165) 



BM^{a) + {I + 5)\1 - MUo')] . 
We can readily demonstrate that Xs < Xc such that physical 
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Figure 5. For versus a of the unaligned case with specified 
5 = 0.2, X = 1-5 S'li'i ^ = the two solid solution curves j/i 
and j/2 are for m = 2 and the one dash-dotted solution curve j/2 
is for m = 1. Here, we use B in the figure to correspond to B in 
the main text. All these solutions increase with increasing a. The 
dashed vertical line represents a = ac- When a becomes greater 
than ac, the branch of yi (m = 2) solution flips and becomes 
negative. In comparison, we see that y2 with m = I is greater 
than J/2 with m = 2. 



solution J/2 > could exist for Xs < X < Xc- However, we 
show presently that Xs > Xc may happen in certain situa- 
tions and the two requirements are incompatible with each 
other; there is thus no physical solution for a positive 1/2 ■ 

The allowed range for x parameter is x > 1- When B = 
1 for the absence of an isopedic magnetic field, the critical Xs 
remains always less than 1 such that the second requirement 
will be satisfied automatically. As the strength of an isopedic 
magnetic field gets stronger, the X> Xs constraint becomes 
indispensable for a positive 1/2 ■ 

Our analysis above bears a strong resemblance to that 
for the aligned case. In fact, should we set the parameter 
Q = for a purely azimuthal propagation of MHD density 
waves, we can approximately obtain nearly the same results 
as those of the aligned case. This naturally suggests that the 
aligned case should be only regarded as a special case of the 
unaligned cases (Lou 2002; Lou & Fan 2002; Lou & Shen 
2003; Lou & Zou 2004). 

Through extensive numerical explorations, we can show 
empirically that both the F-wave solution yi and the S-wave 
solution J/2 increase with increasing a (see Figure 5). From 
the perspective of global stationary MHD density waves, this 
is a physically sensible result. 

4.2.2 Them = l Case 

For m = 1, definition 11561 and recurrence relation 114411 
together give 

1 . .r , . (q2 + 17/4)^/2 



M 



and Ml (a) 



(a2 +5/4) ' ' (a2 +9/4) 

By numerical experiments (see Appendix B) , we have C2 < 
and only j/2 solution being positive. The ratio fi^ /t-i" of j/2 



solution remains always positive with in-phase surface mass 
density perturbations. We also know that j/2 increases with 
decreasing x and increasing B. Numerically, we find that j/2 
solution increases with increasing a. 

The two constraints for sensible j/2 solutions may be 
expressed as 

B 



J/2 > 
corresponding to 

Co > and 



and 



y2 > 



1 



X 



x 



1 + Co > 



respectively (see Appendix B2). 

The mathematical forms of these constraints are too in- 
volved to be shown here. Note that properties of j/2 solution 
with m = 1 are just like those of j/2 in the cases of m > 2 (see 
Fig. 6). As the j/2 solutions of m = 2 and m = 1 are fairly 
similar to each other, the j/2 solution changes smoothly as 
m increases sequentially. 

4.2.3 The m = Case: Marginal Stability 
of Axisymmetric MHD Disturbances 

For m = in definition 115611 and recurrence relation 114411 . 
we have 

2 



M 



(a2 + 1/4) 



and 



Mo{a) = 



(a^ + 9/4) 



(166) 



(a2 + 1/4) (a2 + 17/4)1/2 ' 

On the basis of extensive numerical experiments, we found 
that the five expressions below 



A/'o(c 



1-Moia), 1+M, M+J\fo{a), M+ 



M)(a) 



(1 + 5) ' ' ' ' ' ' ' ' ' (l + S) 

will become negative when a is sufficiently small. As a is 
increased across several different critical values, they will 
in turn become positive. We now proceed to discuss conse- 
quences of this a variation in different ranges separately. 

The analysis is fairly complicated and we summarize 
the basic results here. By imposing the two constraints of 
y > and y > (B/x) — 1, we found that when a < 1.113 
and a > 1.793, there is only one y solution that is physi- 
cally valid under certain conditions. For 1.113 < a < 1.793, 
there is no physically valid solution at all. We then find that 
the valid y solution is characterized by in-phase surface mass 
density perturbations (i.e. a positive f-i^ / jJ." ratio). As before, 
the y solution increases with decreasing x and increasing B 
in both ranges of a < 1.113 and a > 1.793. Numerically, 
we again find that when a < 1.113, the solution of D^ de- 
creases with increasing a. In comparison, when a > 1.793, 
D^ first decreases with increasing a and upon crossing a 
certain point, Dg increases with increasing a. 

For the positive portions of solution without going 
into the extreme, the basic profile of D^ versus a is qual- 
itatively similar to the results of the single SID in Shu et 
al. (2000) and of the composite SIDs in Lou & Shen (2003). 
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Figure 6. The marginal stability curve of D'^ versus a when 
m = 0, B = 2.0, 5 = 0.2 and x = 2.0. Here, we use B in the figure 
to correspond to B in the main text. The solution curve decreases 
with increasing a when a < 1.113. When a > 1.793, the solution 
curve first decreases and then increases with increasing a. There 
exists a minimum for the ring fragmentation curve. 



In fact, our model analysis is more general to include the 

results of Shu ot al. (2000) and of Lou & Shon (2003). 

From the dispersion relation of MHD density waves in 
the WKBJ or tight-winding approximation (Shen & Lou 
2003, 2004b; Lou & Zou 2004), wo sec that the solution 
curve in the versus a profile represents the marginal sta- 
bility curve that separates the regime of stable oscillations 
from the regime of unstable oscillations. As noted by Shu 
et al. (2000), when a < 1.113, the parameter regime under 
the curve in the lower-left corner represents the rotation and 
magnetic field modified Jeans collapse regime. As a stands 
for the radial wavenumber and m stands for the azimuthal 
wavenumber, a small a corresponds to a long-radial wave- 
length MHD density wave disturbance. When a > 1.793, 
the parameter regime above the solution curve represents 
the ring fragmentation regime for large a corresponding to 
short-radial wavelength MHD density wave (see Fig. 7). 

As X parameter increases, the y solution decreases. 
Hence the collapse regime shrinks and the ring fragmentar 
tion regime grows with increasing x (see Fig. 8) . As S param- 
eter increases, the y solution increases. Hence the collapse 
regime grows and the ring fragmentation regime shrinks with 
increasing B (see Fig. 9). 

Physically, x parameter stands for the ratio in the ef- 
fective temperatures of the two SIDs and B represents the 
strength of the isopedic magnetic field. The above results 
mean that when the 'fluid' stellar SID becomes much 'hot- 
ter' than the gaseous isopedic MSID, the stellar SID tends 
to be more stable against collapse for large-scale distur- 
bances but less stable against ring-like fragmentation for 
small-scale disturbances. For a stronger isopedic magnetic 
field, the tendency of collapse becomes easier while the 
ring- fragment at ion becomes more difficult. Here the B pa- 
rameter involves the background equilibrium parameters. 
While we seem to vary only the background magnetic field 



Figure 7. The marginal stability curves of versus a with 
different values of x = 2, 5, 30, respectively, for specified values 
of m = 0, S = 2, and 5 = 1. Here, we use B in the figure to 
correspond to B in the main text. The collapse regime shrinks 
and the ring fragmentation regime grows with increasing x- 
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Figure 8. The marginal stability curves of versus a with 
different values of B = 1, 2, 5, respectively, for specified values 
of m = 0, X = 5, and S = 0.5. Here, we use B in the figure to 
correspond to B in the main text. The collapse regime grows and 
the ring-fragmentation regime shrinks with increasing B. 



strength, the background magneto-rotational equilibrium 
actually changes accordingly. 



5 CONCLUSIONS, SUMMARY AND 
DISCUSSIONS 

For modelling an actual spiral galaxy, our treatment of a 
composite system of MSIDs, while highly idealized still, con- 
tains several physically realistic key ingredients. The formu- 
lation and results of our analysis provide the basic rationale 
for developing useful astrophysical concepts. The two theo- 
rems of Shu & Li (1997) have now been generalized to a com- 
posite MSID system with an isopedic magnetic field. Instead 
of an assumption, the isopedic relation is shown here to be 
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an initial condition maintained by ideal MHD equations. In 
particular, we have taken into account of long-range effects 
of gravity in full without the usual local WKBJ approx- 
imation. We have succeeded in constructing various global 
perturbation configurations, such as bars, spirals, barred spi- 
rals etc. By analogies, we proposed specific transition crite- 
ria among different stationary perturbation configurations. 
With the isopedic magnetic field geometry, we have set the 
stage to further study global patterns of circumnuclear and 
spiral arm galactic winds. We now discuss plausible con- 
sequences and possible applications to spiral galaxies from 
several perspectives. 

On much smaller scales, physical processes involve vari- 
ous star formation activities and their interactions with the 
surrounding ISM. As luminous massive OB star formation 
tends to be more numerous in regions of relatively high gas 
concentration and of stronger magnetic field, the collection 
of numerous young bright stars will then outline a spiral pat- 
tern embedded in relatively high density spiral arms. The- 
oretical studies on the large-scale dynamics of galactic spi- 
ral structure are closely tied to a chain of star formation 
processes and to the physical manifestation of spiral MHD 
density waves. For the global star formation rate in a disc 
galaxy (e.g., Kennicutt 1989; Elmegreen 1994; Silk 1997), 
magnetic field plays an important role in the gaseous disc in 
the context of large-scale instabilities (e.g., Lou & Fan 2000; 
Lou et al. 2001). We would emphasize that should large- 
scale galactic disc instabilities be indeed responsible for the 
global star formation rate, then analyses and simulations 
of a composite disc system involving magnetic field would 
be more realistic and physically meaningful. The oft-used 
axisymmetric stability criterion involving the Q parameter 
is significantly modified in a composite disc system with an 
isopedic magnetic field (Shen & Lou 2003; Lou & Zou 2004). 

Galaxies form and evolve on large scales. To provide 
a dynamic basis for the classification of galactic morpholo- 
gies is one major goal of galactic research (e.g., Berlin et 
al. 1989a, b; Berlin & Lin 1996). The modal approach to 
this problem has gained remarkable progress in explaining 
various configurations of spiral, barred-spiral galaxies with 
observational constraints. For example, a barred galaxy is 
expected to be associated with a relatively heavy disc whose 
halo mass inside the optical radius is small. The potential 
of a rotating bar can give rise to quasi periodic perturba- 
tions in the disc and hence affect the entire configuration of 
a disc. While limited by the stationarity requirement in our 
combined analytical and numerical solution procedure, we 
are able to take into account of long-range effects of gravi- 
tation in full for MHD perturbations in a composite MSID 
system without the deficiency of the usual local WKBJ ap- 
proximation. For example for m = 2, we have constructed 
separately barred configurations (aligned) and logarithmic 
spiral configurations (unaligned). By the principle of linear 
superposition, we might be able to also construct barred log- 
arithmic spiral configurations in a composite MSID system 
by choosing appropriate sets of parameters. In addition to 
theoretically modelling galactic configurations, these solu- 



tions are valuable in testing, initializing and benchmarking 
hydrodynamic and MHD numerical simulation codes for fur- 
ther nonlinear studies. 

As important extension and generalization of the earlier 
theoretical analyses of Shu et al. (2000) on zero-frequency 
(i.e. stationary) aligned and unaligned MHD perturbation 
configurations of an isopedically magnetized SID and of Lou 
& Shen (2003) on stationary aligned and unaligned pertur- 
bation structures in a composite system of two-fluid SIDs, 
we have constructed analytically stationary configurations 
of aligned and unaligned MHD perturbations in a com- 
posite system of a 'fiuid' stellar SID and a 'magnetofluid' 
gaseous MSID. Other closely relevant work include Galli 
et al. (2001), Lou (2002), Lou & Fan (2002), Lou & Shen 
(2003), Shen & Lou (2004a, b), Chakrabarti et al (2003), 
Shen, Liu & Lou (2005), Lou & Zou (2004, 2005). While 
this composite model of one SID and one MSID is highly 
idealized in many aspects, it does contain several necessary 
and more realistic elements that are crucial to understand 
large-scale structures and dynamics of spiral galaxies. 

From the basic fiuid-magnetofluid equations for a com- 
posite system of a stellar SID and a gaseous MSID, we begin 
with an axisymmetric background in a magneto-rotational 
equilibrium. The necessary condition of 1 -f e5 > arises 
here in order to maintain an equilibrium with a positive gas 
surface mass density ('ean. 128^ which in turn puts an up- 
per limit on the magnetic flux under the isopedic condition 
of constant E®/_Bz. One important conclusion of our analy- 
sis is that a constant E®/_Bz is an initial condition sustained 
by ideal MHD equations during disc evolution. Physically, 
the magnetic tension force cannot be too strong to make 
an equilibrium of the gaseous MSID impossible feanl26|l in 
the first place. Massive dark-matter halos tend to stabilize 
disclike galaxies as expected (Ostriker & Feeble 1973; Hohl 
1976; Miller 1978). Although quantitative results for a par- 
tial composite MSID system are not detailed here, we sug- 
gest that the potential of a dark matter halo will certainly 
stabilize a composite MSID system and allow for a stronger 
flux of isopedic magnetic field based on our working expe- 
rience. The B parameter [see definition Ill3ll 1 contains this 
constraint explicitly. 

For MHD perturbations, we derive two sets of linearized 
equations in the stellar SID and the gaseous MSID. Setting 
Lu — 0, we derive the stationary MHD dispersion relation 
with the gravitational coupling. By properly choosing dif- 
ferent potential-density pairs to satisfy the Poisson integral, 
we examine the two classes of aligned and unaligned loga- 
rithmic spiral MHD configurations in parallel. As a result of 
the gravitational coupling between the two SIDs, we obtain 
a quadratic equation in terms oi y = to construct sta- 
tionary MHD perturbation configurations of the composite 
MSID system, namely, equation 110811 for the aligned case 
and equation H155|l for the unaligned logarithmic spiral case. 
The two necessary physical requirements are that both Dg 
and Dg should be positive. 

On the basis of rigorous mathematical derivation and 
analysis, we have reached several conclusions and results. 
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Meanwhile, we have explored parameter regimes numerically 
to reveal a few trends of parameter variations empirically. 

Let us first summarize the basic results of the aligned 
case. For axisymmetric MHD perturbations with m = 0, the 
resulting configuration is simply a rescaling of one axisym- 
metric equilibrium to a neighbouring axisymmetric equilib- 
rium which is certainly allowed but trivial for our purpose. 
For nonaxisymmetric MHD perturbations with m = 1, the 
resulting quadratic equation oi y = is automatically 
satisfied for arbitrary values of y (Shu et al. 2000; Lou & 
Shen 2003). For nonaxisymmetric MHD perturbations with 
m > 2, there are two sets of y = solutions representing 
two different classes of stationary global MHD perturbation 
configurations. For the phase relationships between the sur- 
face mass density perturbations in the two SIDs, the lower 
Ds solutions are featured by an in-phase relation while the 
higher Ds solutions are featured by an out-of-phase rela- 
tion. For a weak isopedic magnetic field, the out-of-phase 
solution always exists (D^ > and Dg > simultaneously), 
while for the existence of the in-phase solution, the condi- 
tion Xs < X < Xc must be met with Xa ™d Xc defined 
by equations 1124^ and H119|l for the aligned case and by 
equations H163|l and H165|l for the unaligned case. As Xs < 1 
may happen as the isopedic magnetic field becomes weak, 
the requirement x > Xa be automatically satisfied by 
the physical constraint x > 1- As the isopedic magnetic flux 
increases across a certain value, the out-of-phase solution 
vanishes while the constraint of the in-phase solution be- 
comes X > Xa with Xs > 1 and with the constraint X < Xc 
being irrelevant. Moreover, we emphasize that in order to 
have Dg > 0, the condition Xa < X rnust be satisfied for 
both aligned and unaligned cases and for both F-wave and 
S-wave solutions; the constraint Xs < X for F-wave solutions 
is satisfied automatically once other requirements are met. 

In the framework of the density wave theory, the case 
of aligned perturbations represents a purely azimuthal prop- 
agation of density wave (Lou 2002) and the stationarity is 
sustained by the advection of a counter SID rotation. In 
general, the class of unaligned logarithmic spiral perturba- 
tions represents MHD density waves with both radial and 
azimuthal propagations where parameter a is an effective 
radial wavenumber. Another special situation is a purely ra- 
dial propagation with m = in the 'unaligned' class. It is 
thus not surprising that with a 7^ in the unaligned class, 
we obtain fairly similar results in parallel with the aligned 
class for m > 2. For example, by setting q = in the un- 
aligned class, we arrive at essentially the same results as in 
the aligned class for m > 2. This justifies our physical inter- 
pretation that the aligned class is merely a special case to 
the unaligned class. For the unaligned class, we obtain two 
solutions of Dg characterized by out-of-phase and in-phase 
surface mass density perturbations, respectively. Other so- 
lution properties of the unaligned class remain qualitatively 
the same as those of the aligned class. 

The unaligned cases with m — 1 and m — need more 
special considerations. For m — 1, only the D^ solution with 
in-phase surface mass density perturbations could be valid 



under certain constraints. For m = 0, only the D^ solution 
with in-phase surface mass density perturbations could still 
exist; when 1.113 < a < 1.793, there is no physical solution 
though. From the local WKBJ dispersion relation, the pro- 
file Ds versus a represents the marginal stability curve (Lou 
& Shen 2003; Shen & Lou 2003, 2004a, b; Lou & Zou 2005; 
Shen, Liu & Lou 2005). In general, when D^ and a are both 
sufficiently small, rotation and magnetic field modified grav- 
itational instabilities for Jeans collapse occur and when D^ 
and Q are both sufficiently large, MHD ring fragmentation 
instabilities appear. Between these two unstable regimes is 
the stable regime of SID rotation for axisymmetric oscil- 
lations. By analytical analyses, we find that all physically 
valid Ds solutions increase with decreasing x and increasing 
isopedic magnetic flux. This is a quite general conclusion 
for both aligned and unaligned cases. The substantial math- 
ematical procedures for these conclusions can be found in 
Appendices A and B. 

Finally, we discuss several possible implications and ap- 
plications of our model analysis. After the theoretical pre- 
diction of the solar wind by Parker (1958, 1963) and the 
subsequent confirmation by spacecraft observations in the 
early 1960s, the more general concept of stellar winds and 
the relevant theories have been developed in parallel (e.g.. 
Burke 1968; Holzer & Axford 1970; Lamers & Cassinelli 
1999). Meanwhile, the concept of galactic winds has grad- 
ually emerged (e.g., Johnson & Axford 1971). Intrigued by 
the almost absence of ISM materials in most elliptical galax- 
ies, Mathews & Baker (1971) also came to the initial concept 
of galactic winds. It is now widely believed that a galactic 
wind (sometimes also referred to as a galactic superwind) 
is powered by starburst activities and supernova explosions 
(e.g.. Chevalier & Clegg 1985; Heckman et al. 1990). Be- 
cause of various observational limitations, the detection of 
a galactic wind has remained a challenge for many years. 
Along with the advance of observational technologies, ev- 
idence for galactic winds have gradually surfaced in the 
multi-wavelength studies. Using the ROSAT data for X- 
ray emissions, Strickland et al. (1997) presented a galac- 
tic wind picture for the nearby galaxy M82. Matthews & 
de Grijs (2004) found the evidence for a galactic wind in 
the edge-on Sbc galaxy UGC 10043 by optical imaging and 
spectroscopy. Melo et al. (2003) reported the detection of 
supergalactic winds in the edge-on starburst galaxy NGC 
4631. Together with galactic-scale outflows, a relevant con- 
cept and phenomenon of superbubbles have been proposed 
and observed (e.g., Tomisaka, Ikeuchi & Habe 1981; Ferriere 
2001). Both phenomena are the results of supernovae, hy- 
pernovae, violent burst activities of massive star formation 
in the densest regions of host galaxies. Galactic winds are 
important processes to channel ISM materials along open 
magnetic flelds into the intergalactic medium, while super- 
bubbles fail to get out of the reign of the host galaxies (e.g., 
Tenorio-Tagle, Silich & Munoz-Tunon 2003). The physical 
reasons why superbubbles cannot become galactic winds re- 
main unclear. 

Based on our model analysis, we here propose that the 
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global structure and geometry of galactic magnetic field 
should be mainly responsible for the coexistence of large- 
scale outflows and supcrbubblcs in disc galaxies. In strong 
analogy to similar phenomena such as global coronal mass 
ejections, violent solar flares from active regions with closed 
magnetic flcld over the solar surface and fast solar wind 
streams coming out of solar coronal hole regions with open 
magnetic flelds, we suspect that bubbles and superbubbles 
may involve closed magnetic fields anchored at the galactic 
disc plane, while large-scale galactic winds pumped by en- 
ergetic sources must be guided by open magnetic fleld lines 
to reach outer regions. 

The model analysis here and those of Lou & Zou (2004, 
2005) and of Shen, Liu & Lou (2005) are complementary 
to each other in terms of overall magnetic field geometries 
associated with discs. In the case of a coplanar magnetic 
field, Parker instabilities inflated by galactic cosmic rays 
and supernova explosions can lead to formation of bubbles 
and superbubbles. In the case of an isopedic magnetic field, 
Velikhov-Chandrasekhar-Balbus-Hawley instabilities, star- 
burst activities and the cosmic-ray pressure together pump 
galactic outfiows from circumnuclear regions and along spi- 
ral arms. We have separately prepared MHD stages in an 
idealized manner for further research on pertinent physical 
processes. In a real disc galaxy, these two different magnetic 
field geometries are intermingled and randomly distributed 
all over the disc. Intuitively, it is natural to imagine that 
outfiows will be stronger at places where open magnetic 
field flux is high. As large-scale dense regions in a galac- 
tic plane, the spiral arms in a galaxy are expected to carry 
more magnetic field flux because of the frozen-in condition 
for magnetic flux in the gas and of small-scale MHD dynamo 
processes. Comparing other regions in the disc (except for 
the nucleus and circumnuclear regions), stronger outflows 
are expected to emerge from spiral axm regions. We refer to 
this scenario as Spiral Arm Galactic Winds. 
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APPENDIX A: 

In the aligned cases of m > 2, we defined two parameters 
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In the quadratic equation of y = Dg for a stationary global 
MHD perturbation configuration 

Civ' + Ciy + Co = , (A3) 
the three coefficients C2, Ci and Co are defined by 
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From quadratic equation IIA3L one can readily derive the 
expression of the determinant A and show that 
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There are thus two real y = solutions for equation l|A3|l . 

For the convenience of analysis, we introduce a useful 
and important Ra parameter, 
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The two following inequalities hold, namely 
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According to the Viete theorem (e.g., John & Horst 1998) 
for a quadratic equation, we know that yi > Ra > 2/2 > — 1 
for C2 > 0, while Ra > y2 > -1 > yi for C2 < 0. Note that 
when C2 < (i.e., X < 0), one must also have Ci < 0. It is 
then clear that yi remains positive and negative for C2 > 
and C2 < 0, respectively. 

To determine the phase relationship between the two 
surface mass density perturbations /i* and /x", we derive 



1 - Him^ 
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where Hi and Gi are separately defined by expressions H98^ 
and ifinnjl . Here, we are primarily interested may solution 
which is at least positive; a negative y solution is unphysical. 
As the case of > always corresponds to yi > Ra > 0, 
it follows that /j.^ / i-i" corresponding to yi solution remains 
always negative for and /j," being out of phase. As 2/2 < Ra 
for any values of C2 , it follows immediately that the ratio of 
/ H" corresponding to j/2 solution remains always positive 
for /I* and ji" being in phase. 

In addition to the physical requirement ol y = > 0, 
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we must also impose Dg > 0. In the purely hydrodynamic 
case of Lou & Shen (2003) without an isopedic magnetic 
field as studied here, Dg > follows from > auto- 
matically. According to the background magneto-rotational 
equilibrium condition of radial force balance, we readily de- 
rive the inequality 



1 > 



Dl > 



which is equivalent to the following inequality 

X 

In other words, besides > 0, we should also require > 
B/x — 1 in the case of a positive RHS. 

We now explore the relations between the two real y so- 
lutions and their variations with specific parameters such as 
X and B. For a chosen parameter, we apply the procedure of 
taking the derivative of the relevant y solution with respect 
to this specific parameter, namely 
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where the prime ' indicates a derivative with respect to this 
chosen parameter. It follows immediately that 
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By simple manipulations, we immediately arrive at 
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where A is the determinant of quadratic equation IIA3II . By 
these expressions, we examine the relations between the two 
y solutions and their dependence on relevant parameters. 

We first discuss solution properties of y\ and 7/2 by vary- 
ing parameter x ^-nd thus have explicitly 
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By expressions IA8II and IIA9L the first derivatives y'l 
and t/2 then take the forms of 
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On the basis of our analysis following equation <A4ll , we 
have y'l < and y'2 < for C2 > 0, and yi > and y'2 < Q 
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for C2 < 0. The basic conclusion is that when yi and j/2 are 
both positive and thus physical, we have both 

^<0 and 5^<0. 
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We next examine the solution properties of yi and 1/2 

as parameter B varies and define explicitly 
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m(m + 2) 



dCo 
dB 



y 



dy 
'dB 



C'2 
Co 



(l + <5) ' 

2m 2(m-|-l) 



(1 + 5) 



+ 



X 

2m(m -I- 1) 



m - 2 



+ 



X 



1 



(1 + 5) 



(1 + 5) 

By rearrangements and manipulations, it is clear that 

Cy^ + Cy + C = 

BiCy'' + Cy + C) + ^22/' +Diy + Do^O , (A12) 
where the three coefficients D2, D\ and Do are defined by 
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Now in the second line of equation IIA12II . we find 
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By the properties of the two real solutions yi and 7/2 and the 
expression above, we can readily show that D2y\ + Diyi -\- 
Do > and -D21/2 + Diy2 -f -Do < 0, respectively. 

With the expression below from the second line of equa- 
tion 1UT2I 
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We note that for a given value of S, parameter S is a function 
of e, while parameter x = Z^/© is yet another function of e 
according to the definition of 0. It is then possible to write 
O as a function of B in the explicit form of 
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together with inequahties derived earUer 
^<0(forC.>0), ^<0 and ^ 
we can clearly demonstrate that 
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^>0 (forC.>0) and ^>0. 

Note that for C2 < and j/i < 0, the case is not considered 
here because of an unphysical negative y — Dl. 

Finally, we proceed to derive analytically the relevant 
constraints of requiring Dl > and Dg > 0, which are 
equivalent to y > and y > (S/x — 1), respectively. 

Let us consider yi first. For C2 > 0, we can readily 
demonstrate that yi > and yi > B/x — ^■ 

We then consider j/2. As we know j/i > B/x — 1 (equiv- 
alent to Dg > 0) and the properties of a quadratic equation, 
we therefore have the following two correspondences 
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For the unaligned cases or the logarithmic spiral cases, we 
first introduce a new parameter A4 as already defined in the 
main text, namely 
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we have the following forms of C2, Ci and Co coefficients 
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(1 + 5)^ ■ ■ X 

The determinant A of quadratic equation IBll can then be 
explicitly shown as 



A = C? - 4C0C2 ^{M + if 



^)] 



"(1 + 5) 



(1 + 5) 



+ 



45AC(q) 



(l + 5)2x 



> 



so that two real y solutions are guaranteed. 

In parallel to what we have done for the aligned case 
in Appendix A, we would like to find out the relationships 
between the y solutions and their variations with relevant 
parameters. Again, we use yi and j/2 to represent the two 
real y solutions defined by 

-Ci±(C?-4CoC2)'/' 



2/1,2 



2C2 



(B2) 



Following the same procedure of Appendix A, we derive 



C'2yl + C(yi +C'a _ C'^yl + C[y^ + C^ 



2/2 



(C? - 4CoC2)i/2 
C'^yl + C[y2 + Ci 



Ai/2 
C'2yl + C[y2 + Co 



(B3) 



(B4) 



(Cf - 4CoC2)i/2 AV2 

where the prime ' denotes the first derivative of a quantity 
with respect to a specific parameter. 

We first discuss the y solution properties with respect 
to the variation of parameter x s.nA define explicitly 



„/ _ dC2 „/ _ dCi 

L^2 = — , '-'1 = — ■ 

dx ox 
It follows immediately that 



C'2 



= 



(X + 1) 



{M + l) 



M 



<^o = — 
OX 



(1+5) 
A/'m(a) 



X' L(l + 5) 

In parallel with the analysis of the aligned case, we here 
introduce a similar parameter Rs in the form of 

' - M+M^{a)/{l + 5) ■ 
It is then straightforward to demonstrate that 
dyi _ {M + l)[M+Nm{a)/{l + 5)] 



dx X2A1/2 

dy2 ^ {M + l)[M+Mm{a)/{l + 5) 

dx X2AV2 



(yi - 



Rs 



y2) 
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We now examine behaviours of y solutions by varying pa- 
rameter B. Following the same procedure, we write 



dCh 
dB 

and derive explicitly 



dCi 
dB 



Co 



1 
g2 



[M+M^{a)] 



M 
M 



dB 



y 



dy 
dB 



(1 + 5) 

" + 



- 1 



(1 + 5) 



It is straightforward to demonstrate that 

dyj_ ^ [M+Mn,{a)][M +M^[a)/{1 + 5)] 
dB B2A1/2 
dy2 ^ [M +N^ {oi)][M+N,^{cx)/{l + 5)] 
dB 



{yi + i)(j/i - 
iy2 + i){Rs 



g2^1/2 

respectively. 

Two important inequalities are summarized below 
(l + MfjVl{a)S 



Rs 



y2) 



CiRt + CiRs + Co 



[M+M5 + Nm{a)Y 



< 



and 



C2{-lf + Ci(-l) + Co = ^^+^^' > . 

X 

The phase relationship between the two surface mass density 
perturbations /i^ and jj." is then given by 

_ l-J/i(m2 + a^ + l/4) 

Rs 



Gi(m2 + a2 + 1/4) 
Afr^ia) + M{1 + 5) [Dl 



Mm{a) (D| + l) ' 

where H\ and d are separately defined by (11481 and 115UII . 
In order to have both Dl and being positive, we still 
have the same two requirements of D^. — y > and — 
y > B/x ~ 1 simultaneously as in the aligned case. 

In the following, we analyze the three situations of m > 
2, m — 1, and m = 0, separately. 

Bl The Cases of m > 2 

Form > 2, the approximation TV™ (q) = l/(m^+Q^ + l/4)^/^ 
is valid (Shu et al. 2000) and we have 
1 



> A/'m(a) > 



1 > X > 



Rs>0 



According to the basic theory of a quadratic equation, we 
infer that yi > i?s > ?y2 > — 1 for C2 > 0, while Rs > y2 > 
— l>yi for C2 < 0. This inference is entirely similar to that 
of the aligned case. And we can get the similar condition for 
j/2 that j/2 > <4> Co > 0. Note also that Rs increases with 
increasing S and a. 

By these analyses, we reach the following conclusions. 
For yi and j/2 being positive and thus physical, we have 



dy2 



< 



dy 



T>0 



dy2 



> 



^<0, 

dx ' dx ' dB ' dB 
The perturbation mass density ratio / ^i" corresponding 
to j/i remains negative with and /i" being out of phase. 



while [I? I ^i" corresponding to 2/2 remains positive with /i® 
and being in phase. The result obtained and the proce- 
dure taken here are just the same as those for the aligned 
case. Moreover, by setting a = in the unaligned case, the 
equation for the unaligned case can be simplified to the same 
equation of the aligned case, implying that the aligned case 
is just a special example of the unaligned case. 

We now indicate consequences of the two physical re- 
quirements 1/ > and y > S/x — 1- By analytical derivations 
for the j/i solution, the two physical requirements 



yi > 



and 



yi > - - 1 

X 



axe simply equivalent to C2 > 0. Meanwhile for the j/2 solu- 
tion, we find that the condition 



is equivalent to 



while the condition 



is equivalent to 

C2f^ 

Vx 



y2 > 



Co >0 



2/2 > - - 1 

X 



+ Ci ( - - 1 I + Co > 



(B5) 



(B6) 



(B7) 



(B8) 



B2 The Case of m 



For m = 1, we have inequalities 

--<M = -(^2 + 5/4) 



< 



and 



0<M(a)=(4±i^<l 



(a2 + 9/4) 

By straightforward algebraic manipulations, we find that 
M + J\fi{a) > and X -I- 1 > 0. As S > 1, C2 remains 
always less than and there must be a negative y solution 
(in fact this y solution must also be less than —1) that cor- 
responds to yi. In this case, only the 2/2 solution may be 
positive in order to be physical. 

Because of j/ij/2 ~ C0/C2 and the solution property of 
yi and C2 derived earlier, the physical requirement of 

y2 > 

is clearly equivalent to 

Co >0 . 

One can readily show the equivalence between the two 
inequalities Rs > and M -|- A/i(a)/(l -I- J) > 0. When 
5 < 0.145307118, Rs remains always positive for arbitrary 
values of a. Furthermore, we know that _Rs > 1/2 > — 1 for 
Rs > 0, while j/2 > -1 > -Rs for _Rs < 0. It then follows that 

fi<0 and ^>0; 
dx dB 

and the ratio of / fi" corresponding to y2 remains always 
positive with /i^ and /i" being in phase. This phase relation- 
ship is just the same as that of y2 solution in the case of 
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m > 2. Here, we have again the following two equivalent 
inequalities between 



2/2 > 
and between 



and 



Co > , 



2/2 > - - 1 
X 



and 



'^'lx~V +Ci(^^-1^ +Co>0 
respectively. 

B3 The Case of m = 

For m = 0, we have 

-8<M = ~ 



Afo{a) = 



(a2 + 1/4) 
(a^ + 9/4) 



<0, 



> 



(a2 + l/4)(a2 + 17/4)1/2 

The situation becomes somewhat involved. First, we identify 

the following five expressions, namely 

1-^, l-Mo{a), 1+M. M+^o{a), M+^'^'^^ 



1 + 5 ' " ' ' ' ' ' 1 + 6 ' 

By extensive numerical computations using the two expres- 
sions of A4 and A/o, we infer several basic results. All these 
expressions are negative when a is small enough. When a 
becomes larger than some critical point, which is different 
for each expression, these expressions become positive. 

For the convenience of analysis, we denote the five criti- 
cal points by ai, a^, as, 04, and as satisfying the conditions 

l-^^^ = 0, l-AAo(a2)=0, 1+M{a3)=0, 

X(a4) +Afo{a4) = , and M{a5) + = , 

respectively. More specifically, numerical calculations indi- 
cate that a2 = 1.113, as = 1.323 and 04 — 1.793. It is found 
that ai < 02 and 04 < 0:5. Actually, ai and 0:5 depend on 
the value of parameter S. For example, when 5 > 3.4, ex- 
pression 1 — A/o(a)/(l + S) remains always positive and ai 
does not exist. In general, our analysis indicates the follow- 
ing ordering for the five Qi parameters, 

ai < a2 < as < a4, < as . 

We now proceed to discuss each a— range separately. 



B3. 1 The Range of < a < ai 

For this a— range to exist, we must require 6 < 3.4. In this 
a— range, we have the following five inequalities 

, A/'o(a) 



(1 + 5) 



<0, l-Afo{a)<0, 1+M<0, 
M+N-o{a)<0, M + ^^<0 



In this case with C2 > and 0.93 > Rs > 0, we have 
2/1 > i?s > 2/2 > — 1. The surface mass density ratio of 



yU^/yti" corresponding to t/i remains positive with /x® and ji" 
being in phase and the ratio of fi^ / fJ," corresponding to 2/2 
remains negative with /x® and /i* being out of phase. We 
further infer 

^<o, ^<0, ^>0, ^>0. 
dx dx ' 9B ' dB 

For the 2/1 solution, we only need to think of the requirement 
2/1 > (B/x) — 1- After some analysis, we can demonstrate 
the equivalence of two inequalities 

2/1 > - - 1 
X 



and 



-If- 



Ci ( - - 1 I + Co < 

.X 



for {B/x) — 1 > Rs- For the 2/2 solution, we found by an- 
alytical derivation that the two inequalities 2/2 > and 
2/2 > (B/x) ^ 1 arc incompatible with each other. Thus in 
this a— range, there is no physical 2/2 solution. 

We note that Rs = for the special case of a = ai , so 
that j/i > _Rs = > y2 > — 1. Here, only the yi solution 
is positive and thus physical with the ratio of jJt^ / being 
positive for in-phase /u® and ^i" . 



B3.2 The Range of a\ < a < a2 
In this a— range, we have inequalities 

1-7Tt4:>0' l-A/'o(a)<0, 1+M<Q, 
(l + o) 

M+Mo{a)<Q, M + ^^(2M<o. 

(1 + <5) 

In this case with C2 > and -0.7444 < 7?^ < 0, we have 
2/1 > J?s > 2/2 > —1, such that 2/2 is negative and unphysical. 

The surface mass density ratio of //-^//i^ of j/i solution is 
positive with /i* and /i^ being in phase. We further infer 



dx 



< 



and 



dyi 
dB 



> . 



Through an analytical analysis, we can demonstrate the 
equivalence of two pairs of inequalities between 

2/1 > - - 1 (B9) 
X 

and 

Ca^^^-l^ +Ci(^^-lj +Co <0 (BIO) 

for {B/x) — 1 > Rs, and between 

2/1 > (Bll) 

and 

Co < . (B12) 

For the special case of a = 02, we still have a positive yi 
solution. However, the conditions yi > and 2/1 > {B/x) " 
1 are incompatible with each other and there is then no 
physical y solution. 
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B3.3 The Range of 02 < a < a-s 

In this a— range, we have the following set of inequalities 
, JVoia) 



and 



(1 + 5) 



>0, 1-A^o(q;)>0, 1 + M<0 
M+Afo{a)<0, M + ^^^<0 



(1+5) 

In this case with C2 > and -1 < < -0.7444, we 
have yi > _Rs > 2/2 > — 1- The 1/2 solution remains negative 
and thus unphysical. The perturbation surface mass density 
ratio of /ti®//u^ for the yi solution is positive with fj,^ and fi^ 
being in phase. We further infer the following inequalities 

fi<0 and ^>0. 
dx dB 

The two inequalities yi > and j/i > (B/x) — 1 are incom- 
patible and there is no physical solution in this a— range. 

When a = as, we have M. = —1 and C2 > 0. It then 
follows that yi — y2 = —1 and there is thus no physical y 
solution for a = 0:3. 

B3.4. The Range of as < a < a4 

In this a— range, we have the following inequalities 

1-7T^>0, l-A/'o(a)>0, 1+M>Q, 
{1 + 6) 

M+Mo{a)<Q, ■^ + -^I^<0- 

With C2 > and Ra < -1, we have -1 > j/i > i?s > 3/2- 
Both j/i and j/2 solutions axe negative and thus unphysical. 

For a = a4, we have C2 = 0, corresponding to the 
divergence point of the quadratic equation; and only one 
negative y solution exists. 

B3.5 The Range of 04 < a < as 

In this a— range, we have the following inequalities 

1-7TT^>0, l-Afo{a)>Q, 1 + M>0, 
(l + o) 

M+Mo{a)>0, >( + ^^°M<o. 

(1 + 5) 

With C2 < and Rs < —1, we have 2/2 > —1 > Ui > Rs- 
The negative j/i solution is unphysical. The perturbation 
surface mass density ratio / of 3/2 solution is positive 
with //® and being in phase. We further infer 

^<0 and ^>0. 

By an analytical analysis, we can demonstrate the equivar 
lence of the following two sets of inequalities between 



and 



^ 1 

2/2 > 1 

X 



C2( ^-1^ +Ci(^^-1) + Co>0, 



and between 



J/2 > 



(B13) 
(B14) 
(B15) 



Co > 



(B16) 



respectively. 

For a = as, wo have C2 < and thus a negative yi 
solution. Only the 2/2 solution has the possibility of being 
positive. The ratio /Li^Z/i" of j/2 solution for a = as is positive 
with //^ and n" being in phase. 

B3. 6 The Range of a > a.j 

In this a— range, we have the following five inequalities 
l-^^^>0, l-A/-o(a)>0, 1+X>0, 

(1 + 5) 

With C2 < and Ra > 0, we have 7?s > 2/2 > -1 > yi- The 
negative yi solution is unphysical. The perturbation surface 
mass density ratio /n" of the 2/2 solution is positive with 
/i* and being in phase. We further infer 

^<0 and ^>0. 

ox oB 

By an analytical analysis, we can demonstrate the equiva- 
lence of the following two sets of inequalities between 



^ 1 

2/2 > 1 

X 



(B17) 



and 

and between 



and 



+ Ci ( ^ - 1 ) + C'o > , (B18) 



2/2 > 



Co > 



(B19) 



(B20) 



respectively. 

Finally, we give a complete summaxy for the m = case 
with different ranges of a values. Physically, we could only 
have one real positive y solution with surface mass density 
perturbations ji^ and being in phase. For < a < a2 = 
1.113 and some constraints, t/i is the only physical solution. 
For 1.113 = a2 < a < a4 = 1.793, no physical solution 
exists. For a > a4 = 1.793 and some constraints, 2/2 is the 
only physical solution. All the physical y solutions increase 
with increasing B and decreasing x- 



